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Abstract 



Dyson-Schwinger equations are integral equations in quantum field theory that describe the Green 
functions of a theory and mirror the recursive decomposition of Feynman diagrams into subdia- 
grams. Taken as recursive equations, the Dyson-Schwinger equations describe perturbative quan- 
tum field theory. However, they also contain non-perturbativc information. 

Using the Hopf algebra of Feynman graphs we will follow a sequence of reductions to convert 
the Dyson-Schwinger equations to the following system of differential equations, 

7[(x) = Pr{x) - sign(s^)7[(x)^ + ^ \sj\-fi{x) xd^jl{x) 

\jeTZ J 

where r G 7?., 7?. is the set of amplitudes of the theory which need renormalization, 7[ is the 
anomalous dimension associated to r, Pr{x) is a modified version of the function for the primitive 
skeletons contributing to r, and x is the coupling constant. 

Next, we approach the new system of differential equations as a system of recursive equations by 
expanding 7i (x) = Z]„>i 7i,„a^"- We obtain the radius of convergence of Z]7i,„a;"/n! in terms of 
that of J2 Pr{n)x'^ /n\. In particular we show that a Lipatov bound for the growth of the primitives 
leads to a Lipatov bound for the whole theory. 

Finally, we make a few observations on the new system considered as differential equations. 
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Chapter 1 



Introduction 



Dyson-Schwinger equations are integral equations in quantum field theory that describe the Green 
functions of a theory and mirror the recursive decomposition of Feynman diagrams into subdia- 
grams. Taken as recursive equations, the Dyson-Schwinger equations describe perturbative quan- 
tum field theory, while as integral equations they also contain non-perturbative information. 

Dyson-Schwinger equations have a number of nice features. Their recursive nature gives them a 
strong combinatorial flavor, they tie Feynman diagrams and the rest of perturbation theory to non- 
perturbative quantum field theory, and on occasion they can be solved, for example However, 
in general they are complicated and difficult to extract information from. 

The goal of the present work is to show how the Dyson-Schwinger equations for a physical 
theory can be transformed into the more manageable system of equations 



where r runs over TZ, the amplitudes which need renormalization in the theory, x is the coupling 
constant, 7i(a;) is the anomalous dimension for r, and Pr{x) is a modified version of the function 
of the primitive skeletons contributing to r, see Chapter [7] for details. 

Chapter [2] discusses the general background with a focus on definitions and examples rather 
than proofs. The approach taken is that Feynman graphs are the primary objects. In an attempt 
to make matters immediately accessible to a wide range of mathematicians and to accentuate the 
combinatorial flavor, the physics itself is mostly glossed over. Readers with a physics background 
may prefer to skip this chapter and refer to existing surveys, such as [T5], for the Hopf algebra of 
Feynman graphs. 

Chapter [3] discusses the more specific background and setup for Dyson-Schwinger equations and 
the insertion operators i?+ on Feynman graphs. Proofs are again primarily left to other sources, 
[l] covers combinatorially similar material for rooted trees. Some important subtleties concerning 

for Feynman diagrams are discussed in more detail in [21j with important results proved in 
[32]. The approach to disentangling the analytic and combinatorial information comes from [24j . 
This chapter leaves us with the following input to the upcoming analysis: combinatorial Dyson- 
Schwinger equations and a Mellin transform for each connected, divergent, primitive graph. The 
former consists of recursive equations at the level of Feynman graphs with the same structure as 
the original analytic Dyson-Schwinger equations. The latter contains all the analytic information. 
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The next four chapters derive (jl.ip expanding upon the discussion in [25] . Chapter H] derives a 
prehminary recursive equation in two different ways, first from the renormalization group equation, 
and second from the Connes-Kreimer scattering-type formula [9j . Chapter [5] reduces to the case 
of single variable Mellin transforms and a single external scale. The Mellin transform variables 
correspond to the different insertion places in the graph, so we refer to this as the single insertion 
place case, though this is only literally true for simple examples. The cost of this reduction is that 
we are forced to consider non-connected primitive elements in the Hopf algebra. Chapter [6] reduces 
to the case where all Mellin transforms are geometric series to first order in the scale parameters 
by exchanging unwanted powers of the Mellin transform variable for a given primitive with lower 
powers of the variable for a primitive with a larger loop number, that is, with a larger number of 
independent cycles. The cost of this reduction is that we lose some control over the residues of the 
primitive graphs. Chapter [7] applies the previous chapters to derive (jl.ip . 

Chapter [8] considers (jl.ip as a system of recursions. It is devoted to the result of [25] where 
we bound the radii of convergence of the Borel transforms of the 7[ in terms of those of Pr- For 
systems with nonnegative coefficients we determine the radius exactly as mm{pr,l/bi}, where pr 
is the radius of the Borel transform of Pr, the instanton radius, and bi the first coefficient of the 
/3-functior0. In particular this means that a Lipatov bouncH for the superficially convergent Green 
functions leads to a Lipatov bound for the superficially divergent Green functions. This generalizes 
and mathematizes similar results obtained in particular cases, such as (/>^, through quite different 
means by constructive field theory [TSj . Both approaches require estimates on the convergent Green 
functions which can also be obtained in some cases from constructive field theory, for example j26] . 

Chapter [9] considers (II. ip as a system of differential equations. We are not able to prove any 
non-trivial results, and so simply discuss some tantalizing features of vector field plots of some 
important examples. More substantial results will appear in 



^This is the physicists' /3-function, see Section [4.11 not the Euler /3 function. 
Lipatov bound for dnul' means that \dn\ < d^n\ for some c. 



Chapter 2 

Background 



2.1 Series 

Definition 2.1. If {an}n>o is a sequence then A(x) = ^n>o'^nX^ is its (ordinary) generating 
function and Yln>o ^^x'^ / nl is its exponential generating function. 

Bold capital letters are used for the ordinary generating function for the sequence denoted by 
the corresponding lower case letters, p will often denote a radius of convergence. 

We will make use of the standard combinatorial notation for extracting coefficients. 

Definition 2.2. If A{x) = J2n>o'^nX^ then [x"]A(x) = an- 

Definition 2.3. Call a power series X]fc>o Gevrey-n if X]fc>o x^a{k) /(fe!)" has nonzero radius 

of convergence. 

For example, a convergent power series is Gevrey-0 and X]fc>o(2^^)'^ is Gevrey-1 due to Stirling's 
formula. Trivially, a series which is Gevrey-n is also Gevrey-m for all m > n. 

Gevrey-1 series are important in perturbative quantum field theory since being Gevrey-1 is 
necessary (but not sufficient) for Borel resummation. Resummation and resurgence are an enormous 
topic which will not be touched further herein; one entry point is [29]. Generally very little is known 
about the growth rates of the series appearing in perturbation theory. They are usually thought 
to be divergent, though this is questioned by some [H], and hoped to be Borel resummable. 

2.2 Feynman graphs as combinatorial objects 

Feynman graphs are graphs, with multiple edges and self loops permitted, made from a specified 
set of edge types, which may include both directed and undirected edges, with a specified set of 
permissible edge types which can meet at any given vertex. Additionally there are so-called external 
edges, weights for calculating the degree of divergence, and there may be additional colorings or 
orderings as necessary. 

There are many possible ways to set up the foundational definitions, each with sufficient power 
to fully capture all aspects of the combinatorial side of Feynman graphs. However it is worth 
picking a setup which is as clean and natural as possible. 
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For the purposes of this thesis graphs are formed out of half edges. This naturally accounts for 
external edges and symmetry factors and permits oriented and unoriented edges to be put on the 
same footing. 

Definition 2.4. A graph consists of a set H of half edges, a set V of vertices, a set of vertex - 
half edge adjacency relations (C T/ x H), and a set of half edge - half edge adjacency relations 
(C if X H), with the requirements that each half edge is adjacent to at most one other half edge 
and to exactly one vertex. 

Graphs are considered up to isomorphism. 

Definition 2.5. Half edges which are not adjacent to another half edge are called external edges. 
Pairs of adjacent half edges are called internal edges. 

Definition 2.6. A half edge labelling of a graph with half edge set H isa bijection H {1, 2, . . . , \H\} 
A graph with a half edge labelling is called a half edge labelled graph. 

2.2.1 Combinatorial physical theories 

Feynman graphs will be graphs with extra information and requirement. In order to define this 
extra structure we need to isolate the combinatorial information that the physical theory, such as 
quantum electrodynamics (QED), scalar (/>^, or quantum chromodynamics (QCD), requires of the 
graph. 

Each edge in the graph corresponds to a particle and a given physical theory describes only 
certain classes of particles, hence the physical theory determines a finite set of permissible edge 
types. For our half edge based setup, an edge type E consists of two, not necessarily distinct, half 
edge types, with the restriction that each half edge type appears in exactly one edge type. An edge 
composed of two adjacent half edges, one of each half edge type in E, is then an edge of type E. 
An edge type made up of the same half edge type twice is called an unoriented edge type. An edge 
type made up of two distinct half edge types is called an oriented edge type. The half edge types 
themselves contain no further structure and thus can be identified with {1, . . . ,n} for appropriate 
n. 

For example in QED there are two edge types, an unoriented edge type, ™" , representing a 
photon, and an oriented edge type, , representing an electron or positroru. At the level of half 
edge types we thus have a half photon, a front half electron, and a back half electron. 

Each vertex in the graph corresponds to an interaction of particles and only certain interactions 
are permitted in a given physical theory, hence the physical theory also determines a set of permis- 
sible vertex types. A vertex type V consists of a multiset of half edge types with 3 < |y| < oo. A 
vertex in a graph which is adjacent to half edges one of each half edge type in V is then a vertex 
of type V. For example in QED there is one type of vertex, . 



The physical theory determines a formal integral expression for each graph by associating a 
factor in the integrand to each edge and vertex according to their type. This map is called the 
Feynman rules, see subsection 12.3.11 On the combinatorial side the only part of the Feynman 
rules we need is the net degree of the integration variables appearing in the factor of the integrand 

^If we chose a way for time to flow through the graph then the edge would represent an electron or positron 
depending on whether it was oriented in the direction of time or not. However part of the beauty of Feynman graphs 
is that both combinatorially and analytically they do not depend on a flow of time. 
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name graph weight 



2 
1 
1 





-1 



Table 2.1: Edge and vertex types in QCD with power counting weights 

associated to each type. Traditionally this degree is taken with a negative sign; specifically for a 
factor N/ D this net degree is deg(Z)) — deg(A^), which we call the power counting weight of this 
vertex type or edge type. 

The other thing needed in order to determine the divergence or convergence of these integrals 
at large values of the integration variables, which will be discussed further in subsection I2.2.4t is 
the dimension of space time. We are not doing anything sophisticated here and this value will be 
a nonnegative integer, 4 for most theories of interest. 

Thus we define. 

Definition 2.7. A combinatorial physical theory T consists of a set of half edge types, a set of edge 
types with associated power counting weights, a set of vertex types with associated power counting 
weights, and a nonnegative integer dimension of space-time. 

More typically the dimension of space-time is not included in the definition of the theory, and 
so one would say a theory T in dimension D to specify what we have called a physical theory. 
Our examples will come from five theories 

Example 2.8. QED describes photons and electrons interacting electromagnetically. As a combi- 
natorial physical theory it has 3 half-edge types, a half-photon, a front half-electron, and a back 
half-electron. This leads to two edge types a photon, ™^ , with weight 2, and an electron, , 
with weight 1. There is only one vertex consisting of one of each half-edge type and with weight 0. 
The dimension of space-time is 4. 

Example 2.9. Quantum chromodynamics (QCD) is the theory of the interactions of quarks and 
gluons. As a combinatorial physical theory it has 5 half-edge types, a half-gluon, a front half- 
fermion, a back half-fermion, a front half-ghost, and a back half-ghost. There are 3 edge types and 
4 vertex types with weights as described in Table 12.11 The dimension of space-time is again 4. 

Example 2.10. t/)^, a scalar field theory, is the arguably the simplest renormalizable quantum 
field theory and is often used as an example in quantum field theory textbooks. As a combinatorial 
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theory it consists of one half-edge type, one edge type, — , with weight 2, one vertex type, X , 
with weight 0, and space-time dimension 4. 

Example 2.11. (j)^ , also a scalar field theory, is another candidate for the simplest renormalizable 
quantum field theory. It is not as physical since to be renormalizable the dimension of space-time 
must be 6, and hence it is not as pedagogically popular. However the Feynman graphs in eft' are 
a little simpler in some respects and so it will be used here in longer examples such as Example 
15.121 </>^ consists of half-edges and edges as in (/>^ but the single vertex type, which has weight 0, is 
3-valent. 

Example 2.12. The final physical theory which we will use for examples is Yukawa theory in 4 
dimensions, which has 3 half-edge types, a half-meson edge, a front half-fermion edge, and a back 
half-fermion edge. The edge types are a meson edge, — , with weight 2 and a fermion edge , 
with weight 1. There is one vertex type, --^ , with weight 0. This example arises for us because of 

2.2.2 Feynman graphs 

Notice that given a graph G, a combinatorial physical theory T, and a map from the half edge of 
G to the half edge types of T, there is at most one induced map from the internal edges of G to 
the edge types of T and at most one induced map from the vertices of G to the vertex types of T. 
Thus we can make the following definition. 

Definition 2.13. A Feynman graph in a combinatorial physical theory T is 

• a graph G, 

• a map from the half edges of G to the half edge types of T which is compatible with the edges 
and vertices of G in the sense that it induces a map from the internal edges of G to the edge 
types of T and induces a map from the vertices of G to the vertex types of T, and 

• a bijection from the external edges of G to {1, . . .n} where n is the number of external edges. 

The final point serves to fix the external edges of G, which is traditional among physicists. 

Lemma 2.14. Let G he a connected Feynman graph with n half edges. Let m he the numher of 
half edge lahelled Feynman graphs (up to isomorphism as lahelled Feynman graphs) giving G upon 
forgetting the lahelling, and let Aut he the automorphism group of G. Then 

m 1 

n! ~ I Aut I 

Proof. Aut acts freely on the n\ half edge labellings of G. The orbits are the m isomorphism classes 
of half edge labellings. The result follows by elementary group theory. □ 

The primary consequence of Lemma 12.141 is that the exponential generating function for half- 
edge labelled graphs is identical to the generating function for Feynman graphs weighted with 
l/|Aut|. 1/1 Aut I is known as the symmetry factor of the graph. Table [22] gives some examples. 

We will be concerned from now on with Feynman graphs which are connected and which remain 
connected upon removal of any one internal edge, a property which physicists call one particle 
irreducihle (IPX) and which combinatorialists call 2- edge connected. Another way to look at this 
definition is that a IPX graph is a unions of cycles and external edges. We'll generally be interested 
in Feynman graphs with each connected component IPX. 
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graph symmetry factor 

bd \ 

1 

2 

Table 2.2: Examples of symmetry factors 

2.2.3 Operations 

For us subgraphs are always full in the sense that all half edges adjacent to a vertex in a subgraph 
must themselves be in the subgraph. 

The most important operations are contraction of subgraphs and insertion of graphs. To set 
these definitions up cleanly we need a preliminary definition. 

Definition 2.15. The set of external edges of a connected Feynman graph is called the external 
leg structure of the Feynman graph. The set of half edge types associated to the external edges of 
a Feynman graph can be identified with at most one edge or vertex type. This edge or vertex type, 
if it exists, is also called the external leg structure. 

Definition 2.16. Let G be a Feynman graph in a theory T, 'y a connected subgraph with external 
leg structure a vertex type V. Then the contraction of 7, denoted G/7 is the Feynman graph in T 
with 

• vertex set the vertex set of G with all vertices of 7 removed and a new vertex v of type V 
added, 

• half edge set the half edge set of G with all half edges corresponding to internal edges of 7 
removed, 

and with adjacencies induced from G along with the adjacency of the external edges of 7 with v. 

Definition 2.17. Let G be a Feynman graph in a theory T, 7 a connected subgraph with external 
leg structure an edge type E. Then the contraction of 7, denoted G/7 is the Feynman graph in T 
with 

• vertex set the vertex set of G with all vertices of 7 removed, 

• half edge set the half edge set of G with all the half edges of 7 removed, 

and with the induced adjacencies from G along with the adjacency of the two half edges adjacent 
to the external edges of 7 if they exist. 

Definition 2.18. Let G be a Feynman graph in a theory T, 7 a not necessarily connected subgraph 
with the external leg structure of each connected component an edge or vertex type in T. Then 
the contraction of 7, also denoted G/7 is the graph resulting from contracting each connected 
component of 7. 
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For example in QED 





Also useful is the operation of inserting a subgraph, which is the opposite of contracting a 
subgraph. 

Definition 2.19. Let G and 7 be Feynman graphs in a theory T with 7 connected. Suppose 7 
has external leg structure a vertex type and let u be a vertex of G of the same type. Let / be a 
bijection from the external edges of 7 to the half edges adjacent to v preserving half edge type. 
Then G o^j ^ is the graph consisting of 

• the vertices of G except for v, disjoint union with the vertices of 7, 

• the half edges of G and those of 7 with the identifications given by /, 
with the induced adjacencies from G and 7. 

Definition 2.20. Let G and 7 be Feynman graphs in a theory T with 7 connected. Suppose 7 

has external leg structure an edge type and let e be an edge of G of the same type. Let / be a 
bijection from the external edges of 7 to the half edges composing e, such that if a is an external 
edge of G then (o, /(a)) is a permissible half edge - half edge adjacency. Then Go^j 7 is the graph 
consisting of 

• the vertices of G disjoint union with the vertices of 7, 

• the half edges of G disjoint union with those of 7, 

with the adjacency of a and /(a) for each external edge a of 7 along with the induced adjacencies 
from G and 7. 

The vertices and edges of G viewed as above are called insertion places. 
For example if 



then there is only one possible insertion place for 7 in G, namely the bottom internal edge e of G, 
and there is only one possible map /. Thus 



then there are 2 possible insertion places for 7 in G, namely the right vertex and the left vertex. 
Let e be the left vertex. Then there are also 4! possibilities for /, however 8 of them give 




On the other hand if 




G Of^g 7 = )00{ 



and 16 of them give 
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Proposition 2.21. 1. Contracting any subgraph ^ of a IPI graph G results in a 1 PI graph. 

2. Inserting a IPI graph 7 into a IPI graph G results in a IPI graph. 

Proof. 1. Without loss of generality suppose 7 is connected. Suppose the result does not hold 
and e is an internal edge in T = G/7 which disconnects T upon removal. Since G is IPX, e 
cannot be an internal edge of G and hence must be the insertion place for 7 in F. However 
then removing either half edge of e from G would disconnect G which is also impossible. 

2. Suppose e is an internal edge in F = G o 7. Removing e removes at least one internal half 
edge of G or of 7 which cannot disconnect either since both are themselves IPX, and hence 
cannot disconnect T. 

□ 

2.2.4 Divergence 

For a IPX Feynman graph G and a physical theory T let w{a) be the power counting weight of a 
where a is an edge or a vertex of G and let D be the dimension of space-time. Then the superficial 
degree of divergence is 

Dl - ^w{e) - ^w{v) 

e V 

where i is the loop number of the graph, that is, the number of independent cycles. Xf the superficial 
degree of divergence of a graph is nonnegative we say the graph is divergent. Xt is the divergent 
graphs and subgraphs which we are primarily interested in. 

The notion of superficial divergence comes from the fact that the Feynman rules associate to a 
graph a formal integral, as will be explained in subsection 12.3.11 the corresponding weights w{a) 
give the degree in the integration variables of the inverse of each factor of the integrand, while 
the loop number £ gives the number of independent integration variables, each running over M.^. 
Thus the superficial degree of divergence encodes how badly the integral associated to the graph 
diverges for large values of the integration variables. The adjective superficial refers to the fact that 
the integral may have different, potentially worse, behavior when some subset of the integration 
variables are large, hence the importance of divergent subgraphs. 

Xn this context we say a theory T (in a given dimension) is renormalizable if graph insertion 
within T does not change the superficial degree of divergence of the graph. 

A theory being renormalizable means more than that the integrals associated to the graphs of 
the theory can be renormalized in the sense of Subsection 12.3.21 Xn fact even if insertion increases 
the superficial degree of divergence, and so the theory is called unrenormalizable, the individual 
graphs can typically still be renormalized. Rather, a theory being renormalizable refers to the 
fact that the theory as a whole can be renormalized, all of its graphs at all loop orders, without 
introducing more than finitely many new parameters. Combinatorially this translates into the 
fact that there are finitely many families of divergent graphs, typically indexed by external leg 
structures. Xn the unrenormalizable case by contrast there are infinitely many families of divergent 
graphs and, correspondingly, to renormalize the whole theory would require infinitely many new 
parameters. 

The interplay of renormalizability and dimension explains our choices for the dimension of space- 
time in our examples. Xn particular 0^, QED, and QCD are all renormalizable in 4 dimensions and 
(f>^ is renormalizable in 6 dimensions. 
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By viewing a divergent graph in terms of its divergent subgraphs we see a structural self- 
similarity. This insight leads to the recursive equations which are the primary object of interest in 
this thesis. 

Another useful definition is 

Definition 2.22. Suppose G is a Feynman graph and 7 and r are divergent subgraphs. Then 7 
and r are overlapping if they have internal edges or vertices in common, but neither contains the 
other. 

2.2.5 The Hopf algebra of Feynman graphs 

The algebra structure on divergent IPI Feynman graphs in a given theory is reasonably simple. 

Definition 2.23. Let TC be the vector space formed by the Q span of disjoint unions of divergent 
IPI Feynman graphs including the empty graph denoted I. 

Proposition 2.24. H has an algebra structure where multiplication m : H <Si H ^ "H is given by 
disjoint union and the unit by I. 

Proof. This multiplication can immediately be checked to be commutative and associative with 
unit I, and to be a linear map. □ 

Another way to look at this is that as an algebra Ti. is the polynomial algebra over Q in 
divergent IPX Feynman graphs with the multiplication viewed as disjoint union. Note that we are 
only considering one graph with no cycles (the empty graph I); from the physical perspective this 
means we are normalizing all the tree-level graphs to 1. 

We will use the notation e : Q ^ H for the unit map e{q) = ql. Also useful is the notation 
Hlin C TC for the Q span of connected nonempty Feynman graphs in Ti. and P\ui : TC TC\in for 
the corresponding projection. That is TCnn is the parts of degree 1. Note that TC is graded by the 
number of independent cycles in the graph, which is known as the loop number of the graph. This 
grading, not the degree as a monomial, is the more relevant in most circumstances. 

The coalgebra structure encodes, as is common for combinatorial Hopf algebras, how the objects 
decompose into subobjects. 

Definition 2.25. The coproduct A : TC TC 'Si TC is defined on a connected Feynman graph F by 

A(r)= 7®r/7 

7cr 

7 product of divergent 
IPI subgraphs 

and extended to 7^ as an algebra homomorphism. 

Note that the sum in the definition of A includes the cases 7 = 1 and 7 = F, since F is divergent 
and IPX, hence includes the terms IiXfr + F^I. Note also that 7 may be a product, that is a disjoint 
union. This is typically intended in presentations of this Hopf algebra, but not always clear. 

Definition 2.26. Let rj : TC ^ Q he the algebra homomorphism with r]{I) = 1 and r]{G) = for G 
a non-empty connected Feynman graph. 

Proposition 2.27. TC has a coalgebra structure with coproduct A and counit rj as above. 
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Proof. We will verify only coassociativity. Calculate (id® A)Ar = J^y t'® A(r/7') = 7'® 
7/7' r/7 where 7' C 7 C F with each connected component of 7' and 7/7' IPI divergent. This 
calculation holds because every subgraph of r/7' is uniquely of the form 7/7' for some 7' C 7 C r. 
Further by Proposition 12.21] and renormalizability each connected component of 7 is IPI divergent, 
so we can switch the order of summation to see that the above sum is simply (A ® id)Ar giving 
coassociativity. □ 

From now on we will only be concerned with the sort of Feynman graphs which appear in TC, 
that is, Feynman graphs with connected components which are divergent and IPI. 

TC is graded by the loop number, that is the first Betti number. 7i is commutative but not in 
general cocommutative. For example in (j)^ theory 

A (^^00^ ) = ^00^ ® I + I ^00^ + 2-Qk + -Q i» • 

Definition 2.28. For fi, f2 ■ T~l ^ H define the convolution /i * /2 = rn{fi /2)A 

We will use the notation id for the identity map TC ^ TC. 
Proposition 2.29. With antipode S : TC ^ TC defined recursively by 5(1) = I and 

5(r) = -F- Yl ^(7) r/7 

7cr 

7 product of divergent 
IPI subgraphs 

on connected graphs, and extended to all of TC as an antihomomorphism, TC is a Hopf algebra 

Proof. The defining property of the antipode is erj = S-k'id = \d-kS. The first equality gives exactly 
the proposition in view of the definitions of A and the second equality is then standard since TC 
is commutative, see for instance [301 Proposition 4.0.1]. □ 

Note that since TC is commutative S is in fact a homomorphism. S is not, however, an interesting 
antipode from the quantum groups perspective since TC is commutative and thus S o S = Id (see 
again [30l Proposition 4.0.1]). 

Definition 2.30. An element 7 of W is primitive if A(7) = ^ + 

A single Feynman graph is primitive iff it has no divergent subgraphs. However appropriate 
sums of nonprimitive graphs are also primitive. For example 

a()oo{ -2X) = ()co{ -2)x( ) 01 + 10 ()oo( 

+^XX ^XX -'XX ^xx 

= ()CO{ -2)x( )0I + I0(^ -^X) 
This phenomenon will be important in Chapter [5l 
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We will make sparing but important use of the Hochschild cohomology of TC. To define the 
Hochschild cohomology we will follow the presentation of Bergbauer and Kreimer [1]. The n- 
cochains are linear maps L : TC ^ TC'^'^. The coboundary operator b is defined by 

n 

bL = (id L)A + ^{-IfAiL + I 

i=l 

where Aj = id • • • id (8> A (8i id (8) • • • id with the A appearing in the ith slot, b^ = since 
A is coassociative and so we get a cochain complex and hence cohomology. The only part of the 
Hochschild cohomology which will be needed below are the 1-cocycles L : TC ^ TC, whose defining 
property bL = gives 

AL= (id(g)L)A + LOl. (2.1) 



2.3 Feynman graphs as physical objects 
2.3.1 Feynman rules 

The information in the Feynman rules is the additional piece of analytic information contained in 
a physical theory, so for us we can define a physical theory to be a combinatorial physical theory 
along with Feynman rules. In the following definition we will use the term tensor expression for a 
tensor written in terms of the standard basis for M.^ where D is the dimension of space-time. Such 
expressions will be intended to be multiplied and then interpreted with Einstein summation. An 
example of a tensor expression in indices /i and is 

P 

where g is the Euclidean metric, /c S M^, /c^ is the standard dot product of k with itself, and ^ is a 
variable called the gauge. Such a tensor expression is meant to be a factor of a larger expression 
like 

where the 7^ are the Dirac gamma matrices, ft is the Feynman slash notation, namely ^ = 'j^k^, 
and m is a variable for the mass. In this example (j2.2p is the integrand for the Feynman integral 
for the graph 




Definition 2.31. Let T be a combinatorial physical theory with dimension of space-time D. Let 
^ be a real variable. Feynman rules consist of 3 maps 

1. the first takes a half edge type (viewed as an external edge), an vector (the momentum), 
and a tensor index fi to a tensor expression in ^, 
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2. the second takes an edge type e, an M-^ vector (the momentum), and tensor indices /x, v for 
each half edge type making up e to a tensor expression in ^, 

3. the third takes a vertex type v and one tensor index //i, /i2, • • • for each half edge type making 
up f to a tensor expression in fii, fj,2, ■ ■ ■■ 

In each case the tensor expressions may depend on 

If there is a non-trivial dependence on ^ in the Feynman rules then we say we are working in 
a gauge theory. QED and QCD are gauge theories. If the Feynman rules are independent of the 
tensor indices then we say we are working in a scalar field theory, (j)'^ and (f)^ are scalar field theories. 
Note that unoriented edges have no way to distinguish their two tensor indices and hence must 
be independent of them. For us the Feynman rules do not include a dependence on a coupling 
constant x since we wish to use x at the level of Feynman graphs as an indeterminate in which to 
write power series. This setup ultimately coincides with the more typical situation because there 
the dependence of the Feynman rules on x is contrived so that it ultimately counts the loop number 
of the graph and so functions as a counting variable. 

Using the Feynman rules we can associate to each graph 7 in a theory T a formal integral, that 
is, an integrand and a space to integrate over but with no assurances that the resulting integral is 
convergent. We will denote the integrand by Int^ and take it over a Euclidean space R-^I^tI where 
D is the dimension of space time and is a finite index set corresponding to the set of integration 
variables appearing in Int^. Then the formal integral is given by 

! Int^ TT d^k 

where D is the dimension of space-time in T and where Int^ and are defined below. 

Associate to each half edge of 7 a tensor index. Associate to each internal and external edge of 
7 a variable (the momentum, with values in M^) and an orientation of the edge with the restriction 
that for each vertex v the sum of the momenta of edges entering v equals the sum of the momenta 
of edges exiting v. Consequently the M- vector space of the edge variables has dimension the loop 
number of the graph. Let be a basis of this vector space. Let Int^ be the product of the 
Feynman rules applied to the type of each external edge, internal edge, and vertex of 7, along with 
the assigned tensor indices and the edge variables as the momenta. 

Note that Int^ depends on the momenta gi, . . . , g„ for the external edges and that these vari- 
ables are not "integrated out" in the formal integral. Consequently we may also use the notation 
Int^(gi, . . . , qn) to show this dependence. The factors associated to internal edges are called prop- 
agators. 

In practice the integrals we obtain in this way are not arbitrarily bad in their divergence. In 
fact for arbitrary A < 00 each will converge when integrated over a box consisting of all parameters 
running from —A to A. 

For example consider 0^ with Euclidean Feynman rules, see [18^ p. 268]. The Feynman rules in 
this case say that an edge labelled with momentum k is associated to the factor l/(fc^ -f-wi^), where 
the square of a vector means the usual dot product with itself and m is the mass of the particle. 
The Feynman rules say that the vertex is associated to —1 (if the coupling constant A was included 
in the Feynman rules the vertex would be associ^ffel^With —A.) Consider 

7 = 

20^ 




oriented from left to right with the momenta associated to the two right hand external edges 
summing to p and hence the momenta associated to the two left hand external edges also summing 
to p. Then the integral associated to 7 is 




1 

{k^ + m?){{p + kY + fn^) 



where df^k = dkQdkidk2dk^ with k = [ko, ki,k2, k^) and squares stand for the standard dot product. 

The above discussion of Feynman rules is likely to appear either unmotivated or glib depending 
on one's background, particularly the rather crass gloss of gauge theories, so it is worth briefly 
mentioning a few important words of context. 

More typically a physical theory might be defined by its Lagrangian £. For example for <j)'^ 

C = \d^ct>d,<P - ^mV' - ^<t>'. 

There is one term for each vertex and edge of the theory and for massive particles an additional 
term. In this case ^d^^(j)dn(j) is the term for the edge of (j)^, —^m^cj)^ is the mass term, and —^(l>^ 
is the term for the vertex. One of the many important properties of the Lagrangian is that it is 
Lorentz invariant. 

The Feynman rules can be derived from the Lagrangian in a variety of ways to suit different 
tastes, for instance directly [HI p. 16], or by expanding the path integral in the coupling constant. 

Gauge theories are a bit more complicated since they are defined on a fibre bundle over space- 
time rather than directly on space-time. The structure group of the fibre bundle is called the gauge 
group. A gauge field (for example the photon in QED or the gluon in QCD) is a connection. A 
gauge is a local section. Choosing a gauge brings us back to something closer to the above situation. 

There are many ways to choose a gauge each with different advantages and disadvantages. For 
the present purpose we're interested in a 1-parameter family of Lorentz covariant gauges called the 

gauges. The parameter for the family is denoted ^ and is the ^ which we have called the gauge 
in the above. The gauges can be put into the Lagrangian in the sense that in these gauges we 
can write a Lagrangian for the theory which depends on ^. For example, for QED in the i?g gauges 
we have (see for example P, p. 504]) 

C = --Ad^A, - d.A^f - ^Ad^A^^f + V^(i7^(5^ - ieA^) - m)^ 

4 Zi^ 

where the are the Dirac gamma matrices. Whence ^ also appears in the Feynman rules, giving 
the definition of gauge theory used above. 

Another perspective, perhaps clearer to many mathematicians is Polyak [27j . 

2.3.2 Renormalization 
Definition 2.32. Let 




be a formal integral. I is logarithmically divergent if the net degree (that is the degree of the 
numerator minus the degree of the denominator) of the integration variables in Int is — Dlfl. I 
diverges like an nth power (or, is linearly divergent, quadratically divergent, etc.) if the net degree 
of the integration variables in Int is — | + n. 
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Let (p be the Feynman rules viewed as map which associates formal integrals to elements of Tl. 
Next we need a method (called renormalization) which can convert the formal integrals for primitive 
graphs into convergent integrals. There are many possible choices; commonly first a regularization 
scheme is chosen to introduce one or more additional variables which convert the formal integrals to 
meromorphic expressions with a pole at the original point. For instance one may raise propagators 
to non-integer powers (analytic regularization) or take the dimension of space-time to be complex 
(dimensional regularization, see for instance [Sj on setting up the appropriate definitions). Then a 
map such as minimal subtraction is chosen to remove the pole part. 

We will take a slightly different approach. First we will set 







(2.3) 



for all r. This is the result which is obtained, for instance, from dimensional regularization and 
from analytic regularization, but simply taking it as true allows us to remain agnostic about the 
choice of regularization scheme. To see the origin of this peculiar identity consider the following 
computation with q E and the square of an element of M.^ denoting its dot product with itself. 



/ 



1 



„r-l 



dx 



il-x[ 



r(r)r(s) Jq (xA;2 + (1 _ x)(/c + q) 



2\r+s 



r(r + S 



r(r)r(s 

T{r + s 



r(r)r(s 
r(r + s 



r(r)r(s 
r(r + s 



r(r)r(s 
r(r + s 



r(r)r(s 

T{r + s 



r(r)r(s 



dxx^ (1 — x) 



d^'k- 



{xk^ + {1 - x){k + qy) 



dxx^ (1 — x) 



s-l 



d'^k- 



dxx'' (1 — x) 



d^'k 



((A; + q{l — x)Y + q'^{x — x^))*"^* 
1 



dxx^ (1 — x) 



s-l 



D 

27rT 



r(#) Jo 



dxx 



r-li 



+q^{x-x'^)y+' 
,1,1 

' '(|/fc|2+g2(2;_x2))-+^ 
\k\D-^ 



dQi 



d\k 



27rf T{r + s-§)T{§] 



2T{r + s) 



(|fc|2 +g2(^_^2))r+s 

1 „ 



dxx 2 



(1-x)^ 



-1-r 







TT-T{r + s - f) 



T{r)T{s) 



D_ 
I 2 ■ 



r(f -r)r(f -.) 

V(p-r- s) 



when 2r + 2s > > 0, > 2r > 0, and D > 2s > 0, and where the first equality is by Feynman 
parameters: 



1 



r(a + /?) 



dx- 



1 



for a, /3 > 



a"6/5 r(a)r(/?) Vo iflx + 6(1 - x))°+/3 

and where dJlfc refers to the angular integration over the unit D — 1-sphere in M^. Now consider 
just the final line and suppose s = 0, then since F has simple poles precisely at the nonpositive 
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integers, is never 0, and 



r(x)r(-x) 



-TT 



X sin(7rx) 

we see that for D > the result is for s = and r not a half- integer. If we view the original 
integral as a function of complex variables r and s for fixed integer D (analytic regularization) , or 
as a function of complex D (dimensional regularization), then by analytic continuation the above 
calculations gives (|2.3p . 

Returning to the question of renormalization, in view of ()2.3p we need only consider logarith- 
mically divergent integrals since by subtracting off in the form of a power of A;^ which is equally 
divergent to the original integral the whole expression becomes less divergent. Logarithmically 
divergent integrals with no subdivergences can then be made finite simply by subtracting the same 
formal integral evaluated at fixed external momenta. 

Let R be the map which given a formal integral returns the formal integral evaluated at the 
subtraction point. In our case then R has as domain and range the algebra of formal integrals 
where relations are generated by evaluating convergent integrals and ()2.3|) . Let </> be the Feynman 
rules, the algebra homomorphism which given a graph G returns the formal integral (piG)- We 
suppose (pil) = 1 and -R(l) = 1. 

If instead we had chosen to use a regulator and corresponding renormalization scheme then (p 
would give the regularized integral of a graph, and R would implement the scheme itself. One such 
example would be dimensional regularization with the minimal subtraction scheme. In that case (j) 
would take values in the space of Laurent series in the small parameter e and R would take such 
a Laurent series and return only the part with negative degree in e. That is R(j){T) is the singular 
part of (pi^): the part one wishes to ignore. Note that in this case -R(l) = 0. The key requirement 
in general is that i? be a Rota-Baxter operator see [H], |15) . 

To deal with graphs containing subdivergences, define Sf^ recursively by Sf^{I) = 1, 

s^r) = -Rim) - E 4(7)i?(0(r/7)) 

7 product of divergent 
IPI subgraphs 

for connected Feynman graphs F extended to all of 7i as an algebra homomorphism. Sf^ can be 
thought of as a twisted antipode; the defining recursion says that Sf^ * R(j) = rj. Use Sf^ to define 
the renormalized Feynman rules by 

4>R = Sf^-k (j). 

When F contains no subdivergences, i;^'_r(F) = (j)(T) — R(j)(T); in view of Subsection 12.3.2] we may 
assume that 0(F) is log divergent and so 0/j(F) is a convergent integral. Inductively one can show 
that (p^ maps TC to convergent integrals. This result is the original purpose of the Hopf algebraic 
approach to renormalization. It gives a consistent algebraic framework to the long-known but ad- 
hoc renormalization procedures of physicists. For more details and more history see for instance 
instance the survey [15] and the references therein. 

These integrals lead to interesting transcendental numbers, but that is very much another story 

m, [22], m- 

Example 2.33. To illustrate the conversion to log divergence and renormalization by subtraction 
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consider the following graph in massless 

k + q 




The Feynman rules associate to it the integral 

1 



<rk- 



J k'^{k + qy 

The factor of l/cp' is there because our conventions have that the graphs with no cycles are all 
normalized to 1. This integral is quadratically divergent and so can not be renormalized by a 
simple subtraction. However we take 



/ 



« 1 



SO 



q^J k^{k + qY q^J {k^ 



2 



q^J {k^Y{k + qY J {k^Y{k + qY 
= -2h-h. 

Each of the two resulting terms are now less divergent. 

To illustrate renormalization by subtraction consider the integral from the second of the above 
terms. As formal integrals (or carrying along the subtraction which we will add below), using the 
same tricks as the calculation earlier this section, 

dS'k^r^ ^ = / <fk C dx- 



(fe2)2(fe + q)2 J - — (^^2 + (1 _ ^-)(^ + q)T 



3 



2 / dxx I d^k 



1 



J {xk^ + {1 - x){k + qff 

1 



{{k + q{l - x)f + q^{x - x"^))"^ 



d^k- 



2 dxx J d^k 

-2 ^ dxx J ^ + ^2(3; _ ^2))3 

"^"'L '''''' L '''*'(|«:P+9*(x-i2))3 
Now consider the result of subtracting at q^ = ji^. By Maple 

\kf \k\' 



rl poo 

I2 - RI2 = 27r^ / dxx / d\k\ 
Jo Jo 



{\k\^ + q\x - (|A;|2 + ^2(a. _ a.2))3 
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= 27r^ J dxx ^-^ log(g^(x - x^)) + ^ log(/x^(x - x^))^ 

= -ylog(^V^^) 

giving us a finite value. 

To finish the example we need to consider the integral 

L=l fd'h ^^-^ 

This integral is linearly divergent so it needs another subtraction of 0. However, this time we only 
need 



/ 



which we can derive from (j2.3p . Write k = k± + /cy where A;|| is the orthogonal projection of k onto 
span(g) and k± is the orthogonal complement, and notice that 



i3 



(4)2 2 y ^{kir 

0-0 = 0. 



So returning to /i, as formal integrals, 



2 / - (^2)2(^ + ^)2 q2 J - (^2)3 



i2 



(fc2)3(A; + g)2 y (fc2)3(fc + g)2- 



The second term is convergent and so needs no further consideration. The first term is now log 
divergent, call it —4/3. Writing k = k^^ + /cy as above, we get 



2 



J {k^)^k + q) 



2 



q^J J -00 H^ + klrikl + ik\\+qy) 

kj\k^\^ 

(fc2 + |fe^|2)3^^2 + (^ , 



00 



2^5/2 /-oo /-oo /j2^4 

' d\k^\ / dfel, " 



r(5/2)7o ' "(A:2 + |fc^|2)3(|fc^|2 + (fc +g)2)' 
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The inner integral Maple can do, and then subtracting at = iJ? the outer integral is again within 
Maple's powers and we finally get a finite answer 

h-Rh = -^7^log{qy^l^). 
lb 

Combining these various terms together we have finally computed / — RI. This completes this 
example. 

Example 2.34. Subtracting in this way also plays nicely with analytic regularization, and is less 
messy on top of it. Consider the example 



'(A;2)i+pi((A; + g)2)i+p2 



Then 



(fc2)l+Pl((fc + g)2)l+P2 J (A;2)l+Pi((fc + g)2)l+P2 (fc2)2+pi+p2 

A;-g((A;2)i+P2 - {(k + qfY+P^^ 



d'^k- 



(fc2)2+pi+p2((^ + g)2)l+p2 



which is merely log divergent and so can be renormalized by subtracting the same integrand at 
q^ = fi'^. This sort of example will be important later on, as we can simply take this integral with 
o = 1 as the Mellin transform which we need in Section [ 



Subtracting off zero in its various forms and subtracting at fixed momenta should not be con- 
fused. The former consists just of adding and subtracting zero and so can be done in whatever 
way is convenient. In the following we will assume that it has been done, and so that all integrals 
are log divergent. The latter, however, we will always explicitly keep track of. It is our choice of 
renormalization scheme and a different choice would give different results. 

2.3.3 Symmetric insertion 

For one of the upcoming reductions we will need to define a symmetric insertion with a single 
external momentum q'^. Let p be a primitive of Tl, not necessarily connected. For the purposes of 
symmetric insertion define the Mellin transform Fp of p (see Section 13. 2p as 

1p1 



where Intp(g2) is the integrand determined by p. We'll renormalize by subtraction at = p? and 
let 

Int-(g2) =Intp(g2) -Intp(//2). 
So define renormalized Feynman rules for this symmetric scheme with subtractions at q^ = /x2 

by 



bR{Bl{X)){q^ / ^?) = I Int;(g2) I ^ |] </>^(X)(-A;2//^') 1 H 



d^h. 
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We have 

<t>R{Bl{X)){q^|^?) = limMX){d^,)FM {{q'/fiY' ' l) , 

where d-p = --§^. 
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Chapter 3 

Dyson-Schwinger equations 



3.1 B 



+ 



For 7 a primitive Feynman graph, Bj_ denotes the operation of insertion into 7. There are, however, 
a few subtleties which we need to address. 

In the closely related Connes-Kreimer Hopf algebra of rooted trees [10] , see Chapter O (F) 
applied to a forest F denotes the operation of constructing a new tree by adding a new root with 
children the roots of each tree from F. For example 



This 1-cocycle property is key to many of the arguments below. The corresponding property which 
is desired of the various appearing in the Hopf algebras of Feynman graphs is that the sum of 
all associated to primitives of the same loop number and the same external leg structures is a 
Hochschild 1-cocycle. 

In the case where all sub divergences are nested rather than overlapping, and where there is only 
one way to make each insertion, a IPI Feynman graph F can be uniquely represented by a rooted 
tree with labels on each vertex corresponding to the associated sub divergence. Call such a tree an 
insertion tree. For example the insertion tree for the graph in Yukawa theory 




B^ in rooted trees is a Hochschild 1-cocycle [10,, Theorem 2] 



AB+ = (id B+)A + B+(g)I. 
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In such cases Bj^ is the same operation as the for rooted trees (with the new root labehed by 
the new graph). So the 1-cocycle identity holds for Bj_ too. 

However in general there are many possible ways to insert one graph into another so the tree 
must also contain the information of which insertion place to use. Also when there are overlapping 
sub divergences different tree structures of insertions can give rise to the same graph. For example 
in (j)"^ the graph 



either at the right vertex or at the left vertex giving two different insertion trees. Provided any 
overlaps are made by multiple copies of the same graph, as in the previous example, then, since 
7 is primitive, the same tensor products of graphs appear on both sides of ()2.ip but potentially 
with different coefficients. Note that this only requires 7 to be primitive, not necessarily connected. 
Fortunately it is possible to make a choice of coefficients in the definition of Bj, which fixes this 
problem. This is discussed in the first and second sections of |2l], and the result is the definition 

Definition 3.1. For 7 a connected Feynman graph define 



where maxf(r) is the number of insertion trees corresponding to F, |X|v is the number of distinct 
graphs obtainable by permuting the external edges of X, bij(7, X, F) is the number of bijections of 
the external edges of X with an insertion place of 7 such that the resulting insertion gives F, and 
(7IX) is the number of insertion places for X in 7. 
Extend Bj_ linearly to all primitives 7. 

Note that -B^(I) = 7. Also with the above definition we have B]^ defined even for nonprimitive 
graphs, but this was merely our approach to make the definition for primitives which are sums; 
now that the definitions are settled we will only consider B]_ for primitives. 

The messy coefficient in the definition of B^ assures that if we sum all B^ running over 7 
primitive IPX with a given external leg structure (that is, over all primitives of the Hopf algebra 




can be obtained by inserting 




into 
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which are single graphs and which have the given external leg structure) , inserting into all insertion 
places of each 7, then each IPI graph with that external leg structure occurs and is weighted by 
its symmetry factor. This property is ^21, Theorem 4] and is illustrated in Example 13. 4[ 

Gauge theories are more general in one way; there may be overlapping sub divergences with 
different external leg structures. Consequently we may be able to form a graph G by inserting one 
graph into another but in the coproduct of G there may be subgraphs and cographs completely 
different from those which we used to form G as in the following example. 

Example 3.2. In QCD 




can be obtained by inserting 




or by inserting 

^ into O. 

This makes it impossible for every B']^ for 7 primitive to be a Hochschild 1-cocycle since there 
may be graphs appearing on the right hand side of (|'2.1|) which do not appear on the left. In these 
cases there are identities between graphs, known as Ward identities for QED and Slavnov- Taylor 
identities for QCD, which guarantee that ^ 5^ is a 1-cocycle where the sum is over all 7 with a 
given loop number and external leg structure. This phenomenon is discussed in |21j and the result 
is proved for QED and QCD by van Suijlekom [32] . 

For our purposes we will consider sets of operators, 

where k is the loop number, r is an index for the external leg structure, and i is an additional index 
running over the primitive graphs under consideration with k loops and external leg structure r. 
In the case where there is only one r under consideration write {i?^'*}*^^Q. Now assume that in this 
more general case, as in QED and QCD, that the required identities form a Hopf ideal so that by 
working in a suitable quotient Hopf algebra we get 

Assumption 3.3. ^*=o -B^'*'*" is a Hochschild 1-cocycle. 

3.2 Dyson- Schwinger equations 

Consider power series in the indeterminate x with coefficients in Ti where x counts the loop number, 
that is the coefficient of x^ lives in the /cth graded piece of 7i. By combinatorial Dyson-Schwinger 
equations we will mean a recursive equation, or system of recursive equations, in such power series 
written in terms of insertion operations . The particular form of combinatorial Dyson-Schwinger 
equation which we will be able to analyze in detail will be discussed further in section 13.31 
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One of the most important examples is the case where the system of equations expresses the 
series of graphs with a given external leg structure in terms of insertion into all connected primitive 
graphs with that external leg structure. More specifically for a given primitive we insert into each 
of its vertices the series for that vertex and for each edge all possible powers of the series for that 
edge, that is, a geometric series in the series for that edge. The system of such equations generates 
all IPI graphs of the theory. 



Example 3.4. For QED the system to generate all divergent IPI graphs in the theory is 

l+2k 

X\ =1 + 



J2 ^'^'^+ 

7 primitive with 
external leg structure 



2k 



J 



X — =l-xB 



X- 



X-^ =I-xB, 



X— X- 



J 



where I7I is the loop number of 7. 

i= the vertex .eri.s. The coofficont of .»,n;. -> is the sum of all IPX QED Feynman 
graphs with external leg structure — ^ and n independent cycles. In QED all graphs have symmetry 
factor 1 so this example hides the fact that in general each graph will appear weighted with its 
symmetry factor. and X~^ are the two edge series. The coefficient of for n > in 

is minus the sum of all IPI QED Feynman graphs with external leg structure and n 

independent cycles. The negative sign appears in the edge series because when we use these series 
we want their inverses; that is, we are interested in the series where the coefficient of consists of 
products of graphs each with a given edge as external leg structure and with total loop number n. 
The arguments to each Bj_ consist of a factor of the vertex series in the numerator for each vertex 
of 7, a factor of the photon edge series in the denominator for each photon edge of 7, and a factor 
of the electron edge series in the denominator for each electron edge of 7. 

To illustrate these features lets work out the first few coefficients of each series. First work out 
the coefficient of x. 



X y =i + xB. 



X- 



-4 

I + xB. ^ (I) + 0(x? 
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X — =I-xB, 



X — X- 



V 



= I - xB^ (I) + 0{x^) 



Next work out the coefficient of x^. 



x-^ 



\ 




= I + XB. 




= I + X 



(l - x-^Q" ) - a:> ^ > ) 




(I) + 0{x^) 




I + a; — C + a; 
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X— =l-xB, 



X- 



I-x- 




+ 0{x^ 



I — x~ 



+ 



X^ =l-xB^ 



X — X- 



I - Xh 



x'^B, 




+ 



+ 0{x^) 



+ 0(x-^ 



The fact that 



and 



appear with coefficient 1 and not 2 is due to the two insertion trees contributing a 2 to the denom- 
inator in Definition 13.11 

By analytic Dyson-Schwinger equations we will mean the result of applying the renormalized 
Feynman rules to combinatorial Dyson-Schwinger equations. These are the Dyson-Schwinger equa- 
tions which a physicist would recognize. The counting variable x becomes the physicists' coupling 
constant (which we will also denote x, but which might be more typically denoted a or depending 
on the theory). The Feynman rules also introduce one or more scale variables Lj which come from 
the external momenta qi and the fixed momentum values //j used to renormalize by subtracting. 
In the case of one scale variable we have L = logg^//i^. See Example 13.51 Note that in the case of 
more than one scale the Lj are not just \ogqf /^'^, but also include other expressions in the qi and 
the Hi, such as ratios of the qi (such ratios are not properly speaking scales, but there is no need 
for a more appropriate name for them since we will quickly move to the case of one scale where 
this problem does not come up). 
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The functions of Lj and x appearing in analytic Dyson-Schwinger equations are called Green 
functions, particularly in the case where the Green functions are the result of applying the renor- 
malized Feynman rules to the series of all graphs with a given external leg structure. 

We can begin to disentangle the analytic and combinatorial information in the following way. 
Suppose we have a combinatorial Dyson-Schwinger equation, potentially a system. Suppose the 
series in Feynman graphs appearing in the Dyson-Schwinger equation are denoted with r £ TZ 
some index set. Denote C the corresponding Green functions. 

For each factor {X'''y in the argument to some Bj, take the formal integrand and multiply it by 
(G^y. For the scale arguments to these C use the momenta of the edges where the graphs of X*" 
are inserted. Then subtract this integral at the fixed external momenta /ij as when renormalizing 
a single Feynman integral. Then the analytic Dyson-Schwinger equation has the same form as the 
combinatorial one but with replacing X'^ and with the expression described above replacing B^. 
Example 13.51 illustrates this procedure. 

In the case with more than one scale the Green functions may depend on ratios of the different 
momenta, and we can progress no further in simplifying the setup. Fortunately, in the case with 
only one scale, which suffices to describe the general case in view of Chapter [5l we can further 
disentangle the analytic and combinatorial information as follows, see |24j for more details. 

Suppose we have a combinatorial Dyson-Schwinger equation and a single scale. For each prim- 
itive graph 7 appearing as a Bj_ we have a formal integral expression 



coming from the unrenormalized Feynman rules. Number the edges, say from 1 to n. Raise the 
factor associated to the ith. edge to 1 + Pi where pi is a new variable. We now have an analytically 
regularized integral which can be evaluated for suitable values of pi . Finally set all external momenta 
to 1. Call the resulting function of pi, . . . , the Mellin transform F^{pi, . . . , pn) associated to 7. 
We are interested in near the origin. 

Then, another way to see the analytic Dyson-Schwinger equation as coming from the combi- 
natorial Dyson-Schwinger equation by replacing X'' with C and Bj^ with F^. The factor with 
exponent pi indicates the argument for the recursive appearance of the X^ which is inserted at the 
insertion place corresponding to edge i. This will be made precise for the cases of interest in the 
following section, and will be motivated by Example 13. 7[ 

Example 3.5. Broadhurst and Kreimer in [5] discuss the Dyson-Schwinger equation for graphs 
from massless Yukawa theory where powers of the one loop fermion self energy > ► are 
inserted into itself. The result is that they consider any graph made of nestings and chainings of 
this one primitive, for example 
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A graph like 



is not allowed. These graphs are in one-to-one correspondence with planar rooted trees. The 
combinatorial Dyson-Schwinger equation is 



X{x) = I - xB^ 



1 



X{x) 



The Mellin transform associated to the single one loop primitive 

k + q 



k 

is, according to the Feynman rules of Yukawa theory, 



1 



k ■ q 



'(A;2)i+Pi((A: + g)2)i+P2 



However we are only inserting in the insertion place corresponding to pi so the Mellin transform 
we're actually interested in is 



1 



F{p) = ^ d^k 



k 



{k'^y+P{k + qf 



£,2 = 1 



Next combine these two facts as described above to get that the Green function satisfies the 
analytic Dyson-Schwinger equation 

k ■ q 



G{x,L) = 1 - — / 



'k'^G{x,log{k^/p^)){k + q)'^ 

where L = log^q^ / p?) and • • • stands for the same integrand evaluated as specified. This is the same 
as what we would have obtained from applying the Feynman rules directly to the combinatorial 
Dyson-Schwinger equation. 

3.3 Setup 

We will restrict our attention to Dyson-Schwinger equations of the following form. 



3.3.1 Single equations 

Fix s G Z. The case s = is not of particular interest since it corresponds to the strictly simpler 
linear situation discussed in [23]. However, to include s = as well, we will make the convention 
that sign(O) = 1. 

Let Q = X^^ . We call Q the combinatorial invariant charge. Applying the Feynman rules to 
Q gives the usual physicists' invariant charge. 
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Consider the Dyson-Schwinger equation 

tk 

X{x) = I-sign{s)Y,Y.^'^+(^Q''^- (3.1) 

fc>l 1=0 

This includes Example 13.51 where s = 2 and there is only one i?+ having k = 1. 

Let Fk^i{pi, . . . , pn) be the Mellin transform associated to the primitive 5^'*(I). In view of 
Chapter [5] we're primarily interested in the case where n = 1 at which point we'll assume that the 
Mellin transforms of the primitives each have a simple pole at p = 0, which is the case in physical 
examples. We expand the Green functions in a series in x and in L (which will in general be merely 
an asymptotic expansion in x) using the following notation 

G(x,L) = l-sign(s)^7fc(a;)L'= jk{x) = ^jk,jX^ (3.2) 

k>l j>k 

The idea is to follow the prescriptions of the previous section to obtain the analytic Dyson- 
Schwinger equation, then simplify the resulting expression by following the following steps. See 
Example 13.71 for a worked example. First, expand G as a series in L. Second, convert the resulting 
logarithms of the integration variables into derivatives via the identity dpy~^\p=o = (— 1)*^ log'^(y). 
The choice of name for the new variable p is not coincidental. Third, switch the order of integration 
and derivation. The result then is a complicated expression in derivatives of the Mellin transforms 
of the primitives. 

However, to avoid the need for additional notation and for appropriate assumptions on the -Ffc^j, 
instead of following this path we will instead define our analytic Dyson-Schwinger equations to be 
the final result of this procedure. 



Definition 3.6. For a single scale p , the analytic Dyson-Schwinger equation associated to ()3.ip 
is 



k 



G{x, L) = 1 - sign(5) ^ ^'^(^' )~"*'"^'^^ ■■■G(.x, )-^isn{^) 



k>l i=0 

{e-^iP^+-+P->^^ -l)F'''\p,,...,p^^ 



Pl = ---=Pnf,=0 



where Uf^ = sign(s)(sA; — 1). 

We only need one subtraction because in view of the discussion at the end of Subsection 12.3.21 
all the integrals of interest are log divergent. 

In view of the following chapters we need not concern ourselves with the complexity of the 
general definition as we will further reduce to the case where there is only one symmetric insertion 
place and a single scale giving 



G{x, L) = l- sign(s) V V x'^Gix, d.pf-'^e-'^P - l)F^^\p) 



k>l i=0 



p=0 



or rewritten 



7-L = Y^ ^'(1 - sign(s)7 • a_p)i-^'=(e-^'' - l)F\p) 



k>l 



(3.3) 

p=0 
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where ^ -U = Y.lkU\ F^{p) = Eto^^'^ip)- 

The connection between the different forms of the analytic Dyson- Schwinger equation and the 
notational messiness of the original presentation can be explained by a motivating example. 

Example 3.7. Let us return to Example 13. 5i The analytic Dyson-Schwinger equation is 



G(x,L) = 1 

where L = \og{q^ / ji^) . 

Substitute in the Ansatz 



X 



k ■ q 



k'^G[x,\og{k'^ / ^i'^)){k + qf 



G{x,L) = l-Y,lk{x)L' 



k>l 



to get 



{k ■ q)ji^ (x) • • • 74 (x) log^(A;V/U^) 



k>l 



^ 7^1 (2;) • • • 7^. (a;) j d'^k 



4 Yl ^h(.^)---iisix) I 



{k.q)log'{k'/fi') 
k'^{k + qf 

{k.q){-lf^%kV^,'rx=o 

k'^{k + qf 



■ diif^'r / d'k 



k- 



' {k'^f+P{k + qf 



p=0 



-1 



2\p 



k>l 



-1 



2\p 



k>l 



.2\p 



Ct'k 



A;- 



d'^kQ 



{BY+P{k + qY 
ko ■ Qo 



p=0 



p=0 



where q = rqo with r € M, r'^ = q'^, q^ = 1 and k = rko 



k>l 



p=0 



using d^y ^\p=Q = (—1)'^ log'^(y). Thus using the notation 7 • f/ = X] IkU^ we can write 

7 • L = x(l - 7 • d-pYHe-'^P - l)F{p)\p=o 

Example 3.8. To see an example of a two variable Mellin transform (a slightly different example 
can be found in p4j) consider again the graph k + p 

7 = 




with the momenta associated to the two right hand external edges summing to p. As an integral 
the Mellin transform of 7 is 



/ 



(A;2)i+Pi((A; + g)2)i+p2 

By the calculations of Subsection 12.3.21 

1 



(A;2)i+Pi((A; + g)2)i+p2 
So the Mellin transform is 



r(i + pi)r(i + p2)' ' r(2-pi-p2) 



F (n n \ - + P2) (2.-p,-p, r(-pi)r(-p2) 

r(l + pi)r(l + P2) T{2-pi-p2) 

Upon subtracting at = ^2 h^qt^ -^g gg^^ 

So the only dependence on q is the dependence on L which is showing up in the correct form for 
Definition 13. 6[ The extra powers of p'^ would get taken care of by the recursive iteration as in 
Example 13.71 

3.3.2 Systems 

Now suppose we have a system of Dyson-Schwinger equations 

X^x) = I - sign(s,) ^ x^B^/''\X''Q^) (3.4) 

k>l i=0 

for r € TZ with TZ a finite set and where 

Q=l[X^{x)~'^. (3.5) 

The fact that the system can be written in terms of the invariant charge Q in this form is typical 
of realistic quantum field theories. For example, in QED (see Example [32 



X 



X-- ] X 



Q 



Suppose a theory T has a single vertex v G TZ with external legs ei G TZ appearing with 
multiplicity rrii, i = 1, . . . ,n where the external legs (made of half-edges types under our definitions) 
are viewed as full edge types, hence as being in TZ, by simply taking the full edge type which contains 
the given half-edge types (hence ignoring whether the edge is the front or back half of an oriented 
edge type). Let val(T;) be the valence of the vertex type v. Then we define 
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For theories with more than one vertex we form such a quotient for each vertex. We are again (see 
section I3.ip saved by the Slavnov- Taylor identities which tell us that these quotients agree, giving 
a unique invariant charge |21l section 2]. Then X^Q^ is exactly what can be inserted into a graph 
with external leg structure r and k loops. 

Proposition 3.9. Suppose Q is as defined as in the previous paragraph. Let G be a IPI Feynman 
graph with external leg structure r and k > loops. Then X^Q^ is exactly what can be inserted 
into G in the sense that we can write X'^Q^ = YYj(z-ji{X^Y^ so that G has tj vertices of type j ft 
j a vertex type and G has —tj edges of type j for j an edge type. 

Proof. In view of (13. for r a vertex type it suffices to prove that we can write 



or 



gfc+(val(r)-2)/2 ^ 

jen 

so that G has tj vertices of type j for j a vertex type and G has —2tj half edges in edge type j 
(including external half edges) for j an edge type. For r an edge type it likewise suffices to prove 
that we can do the same where we define val(r) = 2. 

This holds for some k by viewing a graph as made from a set consisting of vertices each attached 
to their adjacent half edges. 

To see that k is correct note that since G has 1 connected component, e — v + 1 = i, where e 
is the number of internal edges of G, v the number of vertices of G, and i the loop number of G. 
Letting h be the number of half edges (including external half edges) of G we have 

^-. + l-^ = £. (3.7) 
2 2 ^ ^ 

Each Q contributes val(r)/(val(r) — 2) edge insertions and 2/(val(r) — 2) vertex insertions; so each 
Q contributes 

val(r) 2 _ 

val(r) - 2 ~ val(r) - 2 ~ 

to (|3.7p . So if k is so that Q'=+(vai(r)-2)/2 ^Q^j^^^g j^^lf edges and vertices of G as described above 
then 

i = k + '-^^^ + l-'-^ = k. 
So k is the loop number as required. 

□ 

The specific form of Q from (|3.6p will only be used in the renormalization group derivation of 
the first recursion, section 14.11 

Write Fk^i-riPij ■ ■ ■ i Pn^ i-r) ^^^^ the Mellin transform associated to the primitive i3^'*'''(I). Again 
assume a simple pole at the origin. We can then write the analytic Dyson- Schwinger equations as 
in the single equation case. 

Definition 3.10. The analytic Dyson-Schwinger equations associated to (|3.4p are 
G^^(x,Li,...,L,) 
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fc>l i=0 
teTl\{r} 



-L{pi-\ hp. 



sign(st)(s(fc) 
k,i:r 



-sign(st) 



Pl=-=Pnfc^i.^ = 



where the run over the so that the ith factor of is inserted at 



The fohowing notation will be used for expanding the analytic Dyson-Schwinger equations as 
(in general asymptotic) series about the origin, 



G^ix, L) = 1 - sign(.,) E Yk{x)L' Yk{x) = E 

k>l j>k 

In view of the following chapters we will reduce to the case 

G'ix^L) = 

tk 

1 - sign(s^) E E ^''G^'ix, d^pf-'^'' JJ G^{x, d-p)-'j''{e-^P - l)F''^\p) 

je7^x{r} 



(3.8) 



fe>l i=0 



p=0 



or rewritten 



r 

7 • 



L = Y.Y.'^Hl-sig^{sr)Y-d-pf-^'^ n (l-sign(s,)7^-5_, 

jG7^x{r} 



k>l i=0 



(e-^P - l)F'^'*(p) 



(3.9) 



p=0 



where • C/ = X] 7^^^'^- 
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Chapter 4 

The first recursion 



There are two approaches to deriving the first recursion, neither of which is completely self con- 
tained. The first goes directly through the renormalization group equation, and the second through 
the Connes-Kreimer scattering- type formula [9]. 



4.1 From the renormalization group equation 

This is primarily an exercise in converting from usual physics conventions to ours. 

Using the notation of section 13.3.21 the renormalization group equation, see for instance [3 
Section 3.4] or [T7], reads 

4r + P(x)-^ - y Yix) I x(™'(^)~2)/2G'^(x, L) = for I- a vertex type 
oL ox ^-^ } 

e adjacent to / 

d d \ 

h P(x)- 27^(2:) G^(x, L) = for e an edge type 

oL ox J 

where /3(x) is the /? function of the theory, 7^(x) is the anomalous dimension for G^{x,L) (both 
of which will be defined in our notation below), and val(v) is the valence of v. The edge case and 
vertex case can be unified by writing val(e) = 2 and taking the edges adjacent to e to be two 
copies of e itself (one for each half edge making e). Our scale variable L already has a log taken 
so di often appears as in the literature where /i is the scale before taking logarithms. The 
use of a;(™'(^)-2)/2G'^(x,L) in the vertex case in place of what is more typically simply G'"{x,L) 
comes about because by taking the coupling constant to count the loop number rather than having 
the Feynman rules associate a coupling constant factor to each vertex we have divided out by the 
coupling constant factor for one vertex, that is by a;(vai(i^)-2)/2_ ^ result our series begin with a 
constant term even for vertices. 

To see that it makes sense that a vertex v contributes a factor of x(™^(^)~2)/2 recall that that for 
a graph G with one connected component and external leg structure r we have (13. 7p which reads 

2 2 

where h is the number of half edges of G, t the number of vertices, i the loop number, and where 
we take val(r) = 2 for if r is an edge type. Suppose vertices, but not edges, contribute some power 
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of X. Then a vertex v contributes val{v)/2 — 1 to the left hand side of (j3.7p . so it is consistent that 
V also contribute the same power of x. The whole graph G then has x^~^^^^^^^^'^~^ as expected. If 
the power of x associated to a vertex depends only on its valence, then this is the only way to make 
the counting work. 

Returning to /? and 7, define 

(3{x) = dLXcl)R{Q)\L=o (4.1) 

and 

7%x) = -^dLG%x,L)\L=o = llt (4.2) 

for e an edge, that is Se > as discussed in subsection 13.3.21 The factor of x in /3 comes again from 
our normalization of the coupling constant powers to serve to count the loop number (recall from 
the discussion surrounding (j3.6p that Q contributes a 1 to (|3.7p and so, in view of the previous 
paragraph, Q is short one power of x), while the factor of 1/2 in (|4.2p is usual. The sign in (14. 2p 
comes from the fact that our conventions have the Green functions for the edges with a negative 
sign, while the second equality uses the explicit expansion (j3.8p . Note that this /J-function is 
not the Euler f3 function, T{x)T{y)/T{x + y); rather it encodes the flow of the coupling constant 
depending on the energy scale. Another way to look at matters is that the /3-function measures 
the nonlinearity in a theory, specifically it is essentially the coefficient of L in the invariant charge, 
as in the definition above. 

In the case of an edge type e we obtain quickly from (14.11) . (14.21) . (13. Sp . and (13. Sp . that 




s,hiix)x^-Y,ix)]G%x,L). 



In the case of a vertex type v compute as follows. 

0= E 7^(x)jx(™'('')-2)/2G''(x,L) 

\ e adjacent to r / 

= x(™i(^)-2)/2 Ag-(x,L) +x(™'(^)-2)/2/3(x)|-G''(x,L) 
oL/ ox 

+ ^:^Mi:^x(-i-4)/2/3(x)G^(x, L) - x(™'(-)-2)/2 J2 7'(^)G"(x, L) 

e adjacent to r 

oL/ ox 

\ e adjacent to r / 



X— -"- Yl lltix)G^x,L) 



_ ^(val{t;)-2)/2 

2 

e adjacent to r 

from gH), dSSI), and ([33]) 
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Dividing by 3;(™'(^)-2)/2 

using (13.51) and (13.81) we have 

3 en J 
In both cases extracting the coefficient of L'^^^ and rearranging gives 
Theorem 4.1. 



^ ( sign(s^)7[(x) - E 



llix) = ^ I sign(s^)7[(x) - >J |sj|7^(x)xa^ | Yk-iix)- 



Specializing to the single equation case gives 

Theorem 4.2. ^ 

Ik = 77i(x)(sign(s) - |s|x9^)7fc_i(x). 
k 

Note that the signs in the above do not match with [25] because here the sign conventions have 
that the X'^ have their graphs appear with a negative sign precisely if r is an edge type, whereas 
in [25] there was a negative sign in all cases. 

4.2 From 5 

Definition 4.3. Y is the grading operator on TC. Y{^) = I7I7 for ^ TC. 
Definition 4.4. Let 

cTi = dMS*Y)\L=o 
and ^ 

O-n = -:m"-^(f7i (g) • • • (g) f7i)A"-^ 

n! ^ V ' 

n times 

Lemma 4.5. S-kY is zero off Him 

Proof. First S* ★ y (I) = I • = 0. Suppose Li , r2 G W \ QI. Since 5 is a homomorphism and y is a 
derivation, 

S*Y{T,T2) = Y,S{^W,)Y{j'{^'^) 

(E ^(^2)^^(7^') + (E^W)^W')) (E^(^2)7^' 
= 

since by definition 5 * id(ri) = S -k id(r2) = 0. Here we used the Sweedler notation, ^ 7^- ® 7J = 

A(r,). □ 
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Lemma 4.6. 

k 



j=0 
k 

j=0 



where [■] denotes the coefficient operator as in Definition \2.SX 

Proof. The proof follows by induction. Note that both equations read I I = I I when = 0. 
For a given value of A; > the second equality follows from the first for all < ^ < A; using the 
multiplicativity of A and the fact that partitions of k into m parts each part then partitioned into 
two parts are isomorphic with partitions of k into two parts with each part then partitioned into 
m parts. 

Consider then the first equation with /c > 0. By Assumption 13. 3^ for all 1 < i, '^iLo '^+'^ ^ 
Hochschild 1-cocycle. Thus using (j3.4p 



A([x*^]X'-) = A -sign(s,) ^i?f'-([x'=-^]X'-Q^) 

y l<e<k i=0 J 

= -sign(s,) Y X](id»5y'')A([x'=-^]X'-Q^ 



l<i<k i=0 

-sign(.,) Y jZ^Bf^''^[x^~''^X^Q')(^l) 

l<e<k i=0 

i'e / k-e 

= -sign(s,) Y ^(id® 5^"') Y^x^^^'Q^-^ [x^-^-^X'-Q^ 
i<e<k 1=0 \i=o 

+ [x'']X'" I 

= Y(w]^'^Q'~'^-^^Msr) E X^i?f^([x^-^-^]x^-g^) 

j=0 \ l<^<k-j i=0 

+ [x^jX^^ I 

k-l 

= ^[x-^'l^'Q^'-J' ^ [x^'-^jX'' + [x'']X'' I 

j=0 

The result follows. □ 
Note that A(X^) = Efclo^^'Q'' ® (terms of degree k in X"-). 

4.2.1 Single equations 
Proposition 4.7. (T„(X) = sign(s)7„(x) 
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Proof. For n = 1 this appears as equation (25) of [9] and equation (12) of [1]. For n > 1 expand 
the scattering type formula [9l (14)]. The sign is due to our sign conventions, see ()3.2p . □ 

Rephrasing Lemma 14.61 we have 

Corollary 4.8. Suppose Px*^ appears in X with coefficient c and Z Tx^ consists of all terms in 
AX with r on the right hand side. Then Z = cXQ'^ . 

Proposition 4.9. 

(Piin id)AX = X X - sX xd^X 

Proof. The Corollary implies that every graph appearing on the right hand side of AX also appears 
in X and vice versa. Suppose Tx^ appears in X and Z ® Tx^ consists of all terms in AX with F 
on the right hand side. 

By Corollarv l4.8I XQ^ = Z. So in (Piin(8'id)AX we have the corresponding terms P\^^{XQ^)®T . 
Compute 

PU^Q'') = P^,nX + P^r.Q'' 
= PlinX + fcPiinQ 
= PlinX - /csPlinX 

= X - ksX 



Thus 



Theorem 4.10. 



(Piin id)AX = X X - sX xd.^X 



Ik = 77i(x)(sign(s) - \s\xdx)-ik-i{x) 
k 



□ 



Proof. 



7fc = sign(s)(Tj!c(X) by Proposition 14.71 
"°'"^'lm^-i(ai • • • ai)A^-i(X) 



k times 



ai ^^^m^-^( ai0---0ai )A^-^ | A(X) 

\ fc— 1 times 

— ; "i(c7iPiin cjfc_i)A(X) by Lemma 

k 

— sign(s)(Ti(X)cjjt_i(X) — sai{X)xdx(yk^i{X) by Proposition 14.91 

k 

-7i(x)(sign(s)7fc-i(x) - \s\xdx'^k-i{x)) 
k 

□ 
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4.2.2 Systems of equations 
Proposition 4.11. (7„(X'") = sign(sr.)7^(x) 

Proof. The arguments of [9] and [4] do not depend on how or whether the Green functions depend 
on other Green functions, so the same arguments as in the single equation case applied. The sign 
comes from our conventions, see ()3.8p . Note that /3 in [9] is the operator associated to the /5-function 
only in the single equation case, otherwise it is simply the anomalous dimension. □ 

As in the single equation case we can rewrite Lemma 14.61 to get 

Corollary 4.12. Suppose Fx'' appears in with coefficient c and Z Tx^ consists of all terms 
in AX'' with T on the right hand side. Then Z = cX^'Q^. 

Proposition 4.13. 

(Piin «) id) AX' = ® X*" - ^ SjX^ xda.X'' 

3 an 

Proof. As in the single equation case every graph appearing on the right hand side of AX*" also 
appears in X*" and vice versa. Suppose Tx^ appears in X^ and Z Tx^ consists of all terms in AX 
with r on the right hand side. 

By Corollary imiX^Q'^ = Z, and 

= PlinX'- + fcPiinQ 

= PiinX'- -kY, S.jP^^X^ 

jen 

= X^ - A; ^ SjX^ 
jen 

The result follows. □ 



Theorem 4.14. 



Proof. 



^k = -j: { sign(sr)7[(a;)^ - ^ \sj\-fl{x)xdxjl_i{x) 



7^ = sign(sr)o"fc(X'') by Proposition 14.111 
sign(sr) 



k 



-m(cJiPiin o-fc-i)A(X'') as in the single equation case 



sign(sr.)o-i(X'')o-fc_i(X'') - Sjai{X^)xdxSign{sr)(Tk~i{X'') 
by Proposition 14.131 
^ I sign(s^)7[(x)7[_i(x) - ^ \sj\-fl{x)xd^-fl-i{x)) 



□ 
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As in the previous section the signs do not match with [25] because here the sign conventions 
have that the have their graphs appear with a negative sign precisely if r is an edge type, 
whereas in [25] there was a negative sign in all cases. 



4.3 Properties 

The following observation is perhaps obvious to the physicists, but worth noticing 

Lemma 4.15. As a series in x, the lowest term in 7^ is of order at least k. //7i 1 7^ ^ 'l2jen ^jli u 
for ^ = 0, . . . , A; — 1 then the lowest term in 7^ is exactly order k. 

Note that in the single equation case, the condition to get lowest term exactly of order k is 
simply 71,1 ^ 0. 



Proof. Expanding the combinatorial Dyson-Schwinger equation, ()3.ip or (13. 4p . in x we see immedi- 
ately that the x^ term is exactly I. The Feynman rules are independent of x so the x^ term in the 
analytic Dyson-Schwinger equation is 1 — 1 = due to the fact that we renormalize by subtractions. 

Then inductively the 7^ recursion, Theorem 14.11 or 14.141 gives that as a series in x, 7^ has no 
nonzero term before 

y I sign(s^)7[i -{k-l)Y^ \sjhii | 7Li,fc-i 

The result follows. □ 

We also can say that 7^^- is homogeneous in the coefficients of the Mellin transforms in the 
sense indicated below. This will not be used in the following. For simplicity we will only give it in 
the single equation case with one insertion place. 
Expand Fk,i{p) = I]j>-i '^j,k,iP^ ■ 
Recall (1331) 

1-L = Y^ x\l - sign(s)7 • d.pf-'^e-'^f - l)F\p) 

p=0 



7 • L = X x\l - sign(s)7 • d.pf-'^e-'^'^ - l)F\p) 

k>l 

Taking one L derivative and setting L to we get 



71 



X X x\l - sign(s)7 • d.pY-'^pFk,{p)\p=, (4.3) 



k i 



Proposition 4.16. Writing jkj = Yl ^k,j,n,-j^Jj'^h,hM ' ' ' '>^ju,iu,iu for j > k we have that ji + 
\-ju=j-k 

Proof. The proof proceeds by induction. Call ji + ■ ■ ■ + ju the m-degree of 7fcj. First note that 
71,1 = Ei"T-o,i,i from (jMD. 

Assume the result holds for k,j < n. 

Then from Theorem 14.21 or [ 4.101 ^1 „ is a sum over s of terms of the form 

C7£i,ti •••7£„,5„"^, (4.4) 

where £1 + + = s and ti + • • • tu = n — 1. By the induction hypothesis ji^^ti has m-degree 
ti — £i, so ()4.4p has m-degree ^tj — ^£j-|-s = n — 1 — s-|-s = n — 1 as desired. 
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Next from Theorem 14.21 or [ 4.101 o't. j . for k,j < n, is a sum over 1 < -i < n of terms of the form 

C7i,i7fc-i,i-i (4-5) 

By the induction hypothesis 7i^i has m-degree i — 1 and has m-degree j — i — k + 1 so (14. Sp 

has m-degree j — A; as desired. □ 
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Chapter 5 

Reduction to one insertion place 



5.1 Colored insertion trees 

From now on we will need to carry around some additional information with our Feynman graphs. 
Namely we want to keep track of two different kinds of insertion, normal insertion, and a modified 
insertion which inserts symmetrically into all insertion places. Symmetric insertion does not ana- 
lytically create overlapping divergences, but simply marking each subgraph by how it was inserted 
may be ambiguous as in the example below. We will use insertion trees to retain the information 
of how a graph was formed by insertions. 

In examples without overlaps, and even in simple overlapping cases, it suffices to label the 
divergent subgraphs with one of two colors, black for normal insertion and red for symmetric 
insertion. To see that coloring does not suffice in the general case consider the graph 




There are three proper subdivergent graphs; give them the following names for easy reference 




Then if A is red while B and C are black then this could represent A inserted symmetrically into 
XOCX or it could represent B inserted into >(^^ while B itself is made of A symmetrically 
inserted into >{^^ and likewise for C. 

Definition 5.1. A decorated rooted tree is a finite rooted tree (not embedded in the plane) with a 
map from its vertices to a fixed, possibly infinite, set of decorations. 

The polynomial algebra over Q generated by (isomorphism classes) of decorated rooted trees 
forms a Hopf algebra as follows. 
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Definition 5.2. The (Connes-Kreimer) H op f algebra of decorated rooted trees, TCcK, consists of 
the Q span of forests of decorated rooted trees with disjoint union as multiplication, including the 
empty forest I. The coproduct on TCcK is the algebra homomorphism defined on a tree by 

A{T) = Y,PciT)0RciT) 

c 

where the sum runs over ways to cut edges of T so that each path from the root to a leaf is cut at 
most once, Rc{T) is the connected component of the result connected to the original root, and Pc{T) 
is the forest of the remaining components. The antipode is defined recursively from 5* id = ery (as 
in the Feynman graph situation). 

See [To] for more details on TCck- Insertion trees are decorated rooted trees where each element 
in the decoration set consists of an ordered triple of a primitive of TC (potentially a sum), an insertion 
place in the primitive of the parent of the current vertex, and a bijection from the external edges 
of the Feynman graph to the half edges of the insertion place. The second and third elements of 
the triple serve to unambiguously define an insertion as in Subsection 12.2.31 Often the insertion 
information will be left out if it is unambiguous. 

Definition 5.3. For a IPX Feynman graph G in a given theory let F{G) be the forest of insertion 
trees which give G. 

From F{G), or even just one tree of F{G), we can immediately recover G simply by doing 
the specified insertions. The result of the insertion defined by a particular parent and child pair 
of vertices is unambiguous since all the insertion information is included in the decoration. The 
choice of order to do the insertions defined by an insertion tree does not affect the result due to 
the coassociativity of the Feynman graph Hopf algebra. 

Extend F to F : TC ^ TCck as an algebra homomorphism. In fact it is an injective Hopf algebra 
morphism by the following proposition. 

Proposition 5.4. F{A{G)) = A(F(G)). 

Proof. Let 7 be a (not necessarily connected) divergent subgraph of G. Since G can be made by 
inserting 7 into G/'j, then among F{G) we can find each tree of -^(7) grafted into each tree of 
F{G/^). Cutting edges where -^(7) is grafted into F{G/'y) we see that -^(7) F{G/^) appears in 
A(F(G)). The coefficients are the same since each insertion place for 7 in G/^ which gives G we 
have a grafting with this insertion information and vice versa. Finally every cut of F{G) consists 
of a forest of insertion trees, which by doing the insertions gives a divergent subgraph of G. The 
result follows. □ 

Now we wish to extend this situation by coloring the edges of the insertion trees. 

Definition 5.5. Let T be a decorated rooted tree with edge set E. Define an insertion coloring 
map to be a map f : E ^ {black, red}. If T is an insertion tree when call T with / a colored 
insertion tree. 

Definition 5.6. For a colored insertion tree define the coproduct to be as before with the natural 
colorings. 
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To translate back to Feynman graphs think of the edge as coloring the graph defined by the 
insertion tree below it. The result is a Feynman graph with colored proper subgraphs. The 
coproduct in the tree case forgets the color of the cut edges. Correspondingly in the Feynman 
graph case the color of the graphs, but not their subgraphs, on the left hand sides of the tensor 
product are forgotten. 

Proposition 5.7. Colored insertion trees form a Hopf algebra with the above coproduct which 
agrees with TicK upon forgetting the colors. 

Proof. Straightforward. □ 

Call the Hopf algebra of colored insertion trees Tic. In view of the above R : 7i ^ Tic by taking 
R : 7i ^ TLcK and coloring all edges black. 

Analytically, black insertion follows the usual Feynman rules, red insertion follows the symmetric 
insertion rules as defined in subsection 12.3.31 

Definition 5.8. For 7 a primitive element of 7i or Tic, write i?^ : Tic Tic for the operation of 
adding a root decorated with 7 with the edges connecting it colored red. Also write Bj, : Tic Tic 
for ordinary insertion of Feynman graphs translated to insertion trees with new edges colored black. 
Note that this is not the usual i?+ on rooted trees in view of overlapping divergences. 

When working directly with Feynman graphs i?^ corresponds to insertion with the inserted 
graphs colored red and no overlapping divergences. 

Another way of understanding the importance of Definition 13.11 and Theorem 13.31 is that 
Z^iLo B^^'^ is the same whether interpreted as specified above by 5+ on Feynman graphs translated 
to Tic, or directly on Tic simply by adding a new root labelled by 7 and the corresponding insertion 
places without consideration for overlapping divergences. 

Lemma 5.9. i?^ is a Hochschild 1-cocycle for Tic. 

Proof. The standard i?+ of adding a root is a Hochschild 1-cocycle in TicK, see \10\ Theorem 
2]. Edges attached to the root on the right hand side of the tensors are red on both sides of the 
1-cocycle identity. The remaining edge colors must also satisfy the 1-cocycle property which we can 
see by attaching this information to the decoration of the node which is further from the root. □ 

5.2 Dyson- Schwinger equations with one insertion place 

To reduce to one insertion place we need to show that we can write Dyson- Schwinger equations 
in which only involves i?+s but which, results in the same series which contains only black 
insertions. We can achieve this recursively, while viewing Ti ^ Tic. 
Suppose our combinatorial Dyson-Schwinger equation is as in (j3.4p 




k>l i=0 



Then, using [•] to denote the coefficient operator as in Definition 12. 2^ define 
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n-1 



k=l 

n ''k n— 1 



fc=l 4=0 fc = l 

In order to know that the gj^ are weh defined we need to know that they are primitive. 
Proposition 5.10. is primitive for r eTZ and n > 1. 

Proof. First note that is primitive for any and the sum of primitives is primitive, so q[ 

is primitive for each r ^ TZ. 
Then inductively for n > 1 

n '-fc n— 1 



k=l 1=0 k=l 

n '■k n—1 

+ ^Y1 B'l!'''{[x''-'']X'Q'') I - ([x'^-'^lX'-g^ 

k=l i=0 k=l 

n n — k 

A;=l i=0 £=0 
n—1 n—k 



fc=l £=0 



e=i k=i \ \i=o 

n-l 



Theorem 5.11. 

X"- = 1- sign(s^)^x'^i?5(X' Q^). 
fe>i 

Proof. The constant terms of both sides of the equation match and for n > 1 

n 

-sign(s^)[x"]^x'=i?5(X''Q'=) = -sign(s^)^x'=i?J([x"-'=]X''Q'=) 

k>l k=l 

n-l 

= -sign{sr)ql - sign(s,) ^ x'^i^J ([x"-'=]X^Q'=) 



fc=i 
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□ 



□ 

The interpretation of the Theorem is that we can reduce to considering only red insertion, that 
is to a single symmetric insertion place. 

In simple cases we can avoid the not only the insertion trees, but also the subgraph coloring, 
and literally reduce to a single insertion place in the original Hopf algebra. However this cannot 
work with different types of insertions or with vertex insertions where each vertex can not take 
an arbitrary number of inserted graphs. Consequently such simple examples can only arise with a 
single type of edge insertion as in the following example. 

Example 5.12. Suppose we have the Dyson-Schwinger equation 

X = l-xBl^ (i^). 

where we insert into both internal edges. In this case we need not resort to red insertion in order 
to reduce to one insertion place. 
Let 




where we only insert into the bottom edge and let 

X. = l-xB'.(i,) 

Then to order we have that 

X = l-x-^ ~x^-^d^ -x^(^—%^ +1^©^ +1^^^ 

and 

^i=i-4o -^i^ -^'G^ +\^) 

so 1 A 1 1 

02 = and Qz = - -G^ ^oo^ 

^ ^ 8 16 16 

where in the first graph of we insert only in the bottom edge of the bottom inserted bubble, 

in the second graph we insert only in the bottom edge of the leftmost inserted bubble, and in the 
third graph wc insert only in the bottom edge of the rightmost inserted bubble. 
Note that is primitive. Let 

The order x^ we have 

x=i-x-0 -x'l^d^ -^^{l^ + \^^ + \^^^ 

2^2 V8 2 4 y 
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and 



.^4/1^^ +1^^ +1^^ +1^^ +1^^ 
V8 4 2 8 4 

+1^^ ^ 
8 4 4 / 



^2 \2 8 y 

^^^L +1^^^ +^^$^ 
8 4 2 8 8 

8 8 / 



where the first 2 Unes come from inserting X2 into q'l and the third hne comes from inserting X2 
into qs- 

Consequently let 

Qyi = — — — —'000^ — 1 + 1 — 'o®^ + 1 — /^o^ 

^ 8 8 4 8 8 

which we can check is primitive. Continue likewise. 
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Chapter 6 

Reduction to geometric series 



6.1 Single equations 



Let D = sign(s)7 • d-p and Fk{p) = X]i=o ^k,iip) so the Dyson-Schwinger equation (|3.3p reads 



7 • L = ^ x'il - Df-^^e-^" - l)F,{p) 
k>l 



p=0 



Only terms Ux^ with k > j > 1 occur by Lemma 14.151 so this series hes in (M[L])[[a;]]. Then we 
have the fohowing 

Theorem 6.1. There exists unique rk,rk,i GM, /c>l, l<i<fc such that 

j2xH'i--D)'-'\e-'^'-mip) 



p=0 



j:ai-d)'-v^'-i)|-7j^+ e 



P(l - P) , 



<i<k 



p=0 



Proof. For £ > the series in x 



x\l-Df-'^pX=^ 



has no term of degree less than k + 1 since 7i(x) has no term of degree less than i by Lemma 14.151 
It follows that 

x^{l- Df-'^'ie-'^P -I) ^ 



and 



\l- Df-'^{e-^P -I 



P(l-P) 
U 



-Lx^ + 0{x^+^] 



p=0 



P 



V+^x^ + 0{x 



fc+l^ 



p=0 



Now expand 



i=-i 
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and define r„, rn,i recursively in n so 



A; 



p=0 



+ E 



l<i<k 



+ 0{x 



p=0 



This is possible since as noted above the coefficient of x"' in x^{l — D)^ '^^{e — l)Ffc(yo)|p=o 
is a polynomial in L with degree at most n — 1. □ 

The meaning of this theorem is that we can modify the Mellin transforms of the primitives to 
be geometric series at order L. The higher powers of p in the Mellin transform of a primitive at k 
loops become part of the coefficients of primitives at higher loops. Note that there are now terms 
at each loop order even if this was not originally the case. 

Example 6.2. Consider the case s = 2 with a single Bj^ at order 1 as in Example 13.51 Write 

oo 

i=-i 

then computation gives 

Tl =f-l 

r2 =/-i - /-i/o 
r2,i =0 

n =2/!i + /!i(-4/o + /i) + f^ifi 

r3,i = -f-i + f-ik 
r3,2 =0 

r4 =2/li + /!i(-12/o + 6/i - /2) + /!i(9/2 - 3/o/i) - f^if^ 

U,i = - /li + /'i(6/o - 2/i) - 3/2 J2 
7 7 
D 

rs = - 10/5, + /li(-6/o + 18/i - 8/2 + /s) + /!i(40/2 - 32/o/i + 4/0/2 + 2/^) 
+ /!i(-16/3 + 6/2/i)+/-i/o 



These identities are at present still a mystery. Even the coefficients of f^i in do not appear 
in Sloane's encyclopedia of integer sequences [28] in any straightforward manner. In the case 



p(l-p)(2-p)(3-p)' 
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as in the 0^ example from [5], the above speciahzes to the also mysterious sequence 

1 
5 

r2 = - ^ -2,1 =0 

14 -5 
r-s = - ^ -3,1 r3,2 =0 

563 -173 -35 

13030 

-5 



69 

194178 



re 



6 



11 



Note that even if the coefficients of the original Mellin transforms are all of one sign the 
unfortunately not be so. 

6.2 Systems 

As in the single equation case we can reduce to geometric series Mellin transforms at order L. 
Theorem 6.3. There exists unique r^,r^^ G M, A; > 1, l<i<k, j£TZ such that 



tk 

x^{l - sign{sr)Y ■ d-pf-'^-^ W (1 - sign(sj)7^a_p)-^^'=(e--^'' - l)F^'\p) 

k>l i=0 j<zTi-^[r} 

= EE^'(l-sig^^(^-)^'-^-p)'"'^' n (l-sign(s,)7^5_p)-^^*^ 



p=0 



k>l 1=0 jG7^\{r} 



p=Q 



Proof. The proof follows as in the single equation case with the observation that for i >0 



X 



still has lowest term x^~^^. 
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Chapter 7 

The second recursion 



7.1 Single equations 



Reducing to geometric series Mellin transforms at order L allows us to write a tidy recursion for 
7i. Again let D = sign(s)7 • d-p and Fk{p) = Yl^iLQPk,i{p)- By Theorem 16.11 we have 



p=0 



Taking the coefficients of L and gives 
So 



p=0 



p=0 



71 + 272 = ^xHl - Df-'\-ru - 2rfe,i) = Y.p{k)x^ = P{x) 
k P=o k 



where p{k) = — — 2rk^i. Then from Theorem 14.21 or [ 4.101 

71 = P{x) - 272 = P{x) - 7i(sign(s) - \s\xdx)-ii 

giving 

Theorem 7.1. 

71 (x) = P{x) - 7i(x)(sign(s) - \s\xdx)-ii{x) 
or at the level of coefficients 

n-l 

7i,n = p{.n) + ^{\s\j - sign(s))7ij7i,„_j 



(7.1) 
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Notice that in defining the and r^^j we only used a geometric series in the first case. Specif- 
ically, we used l/(p(l — p)) for but 1/p for r^^j. We could have used in all cases; then 
one L derivative would give 71 (a;) = "^r^x^ so all the information of 71 is in the r^, we learn 
nothing recursively. The choice of a geometric series at order L was made to capture the fact that 
conformal invariance tells us that the Mellin transform will be symmetrical when p ^ 1 — P) and 
it also entirely captures examples such as the Yukawa example from [5j and Example 13.51 On the 
other hand choosing to use a geometric series for the j as well would not have resulted in a tidy 
recursion for 71 using these techniques. We hope that the choice here gives an appropriate balance 
between representing the underlying physics and giving tractable results all without putting too 
much of the information into P{x). 

Another important question is how to interpret P{x). In cases like the Yukawa example of [5] 
where the various reductions are unnecessary, then P{x) is simply the renormalized Feynman rules 
applied to the primitives. In that particular example there is only one primitive, and P{x) = cx 
for appropriate c. In the general case we would like to interpret P{x) as a modified version of the 
renormalized Feynman rules applied to the primitives. For the first reduction this is a reasonable 
interpretation since that reduction simply makes new primitives, either within the Hopf algebra of 
Feynman graphs or more generally. For the second reduction the idea is that the geometric series 
part of each Mellin transform is the primary part due to conformal invariance. At order L the rest 
of the Mellin transform gets pushed into higher loop orders, while at order L? the reduction is a bit 
more crass. This information together gives the and the r^^i and hence gives P{x). So again, in 
view of the previous paragraph, we view this as a modified version of the Feynman rules applied to 
the primitives. 

7.2 Systems 

Theorem 16.31 gives us 

fc>i i=o jen-^ir-} 



p=0 



As in the single equation case we can find tidy recursions for the 7[ by comparing the coefficients 
of L and in the above. We get 



l{ = -^x''{l-sign{srh'.d^,)'-'^'' n (l-sign(s,)7^'-a_,)-^^'=^|p=o 



and 



272 = E - sign(s,)7'^ • d^p)'-''-' n (1 - sign(^j)7' " d- 

k iG7e\{r} 

P^k 
1-p 



\—s.ik 
p) 



2rfc,i \p=o 
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k>l 

Thus letting p'^{k) = —r^ — 2r'^ ^ and using the first recursion (Theorem 14. II or 14. 14p 

7[ = ^pHk)x'' - 27^ = ^/(A;)x'= - sign(s,)7[(a;)^ + Yl \sihiix)x9^^lix) 
k>l k>l jell 

giving 

Theorem 7.2. 



7i = Yp''{k)x'' - sign{sr)-fi{xf + ^ \sj\-fl{x)xdxjl{ 
k>i jen 

or at the level of coefficients 

n—l n—1 

li,n = + Yi\sr\i - sign(s^))7[_i7[^„_i + ^ ^ilsjlihi^n-nli 



1=1 j&n 1=1 



7.3 Variants 



The value of the reduction to geometric series is that ifF(p) = r/(p(l— p)) then p^F(p) = pF{p) — l. 
However this reduction is rather crass, particularly for higher orders of L, so it is worth considering 
other special forms of F as in the following example. 

Example 7.3. Consider again the 0^ example from [5] as setup in Example 16. 2[ We have s = 2 
and 



p{l - p){2 - p){3 - pY 



so 



pF{p) ^ 



(l-/>)(2-p)(3-p) 



ay (i-p)(2-p)(3-p)y 

l + p^Fip)-^p'Fip) + lp'Fip) 
6 6 6 



This gives that 



71 =-x(l-7-a_p) ^pF{p)\p=o 

=^x{l - 7 • 5_p)-il|,=o - x(l - 7 • d-,)-'p^Fip)\p=o 

+ - 7 • d.p)-'p^Fip)\p=o - ^xil - 7 • a-p)-V'i^(p)|p=o 

=7- - 272 - II73 - 474. 
6 
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In view of Theorem 14.21 or 14.101 which in this case reads 

Ik = -^71 (2;) (1 - 2x9^)7fc_i(x), 

we thus get a fourth order differential equation for 71 which contains no infinite series and for which 
we completely understand the signs of the coefficients. 
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Chapter 8 

The radius of convergence 



8.1 Single equations 

We see from the second recursion, Theorem 17.11 that if '^p{k)x^ is Gevrey-n but not Gevrey-m 
for any m < n, then 71 is at best Gevrey-n. 

Of most interest for quantum field theory apphcations is the case where only finitely many p{k) 
are nonzero but all are nonnegative and the case where p{k) = c^kl giving the Lipatov bound. In 
both cases Y^p{k)x^ is Gevrey-1. Also for positivity reasons we are interested in s > 1 or s < 0. 
Thus for the remainder of this section the following assumptions are in effect. 

Assumption 8.1. Assume |s| > 1 or s < 0. Assume p{k) > for k > 1 and 

fc>i 

has radius of convergence < p < cxd and is not identically zero. 

Under these assumptions 71 is also Gevrey-1 and the radius is the minimum of p and l/(sai) 
(where we view l/(sai) as +00 in the case oi = 0) the proof of which is the content of this section. 

Definition 8.2. Let = ji^n/nl, A{x) = J2n>i (^nX^, and let pa be the radius of convergence of 
A(x). 

Then ai = 71^1 = p{l) and 



.?^+l:(H.-.g„(.))(" 

j=l v-^ 

, ^ , 71-1 

p[n 



3 



1 "~ /n\ 
+ l^'^{\s\j - s\gii{s) + \s\{n - j) - sign{s)){ \ ajan- 

.7=1 ^"'^ 

n-l / s 



P^""^ + f\s\!l - sign{s))y2 i'') a.an^j (8.1) 



n! V '2 

Inductively, we see that ai , 02 , ... are all nonnegative 
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Note that if s = 1, p{l) > 0, and p{n) = for n > 1 then oi = p{l), ffln = for n > 1 solves the 
recursfon. In this case pa = p = oo, but < l/(|s|ai) < oo. This boundary case is the only case 
with this behavior as we see in the following Proposition. 

Proposition 8.3. Suppose that either s ^ 1, orp[n) > for some n > 1. Then pa < min{p, l/(|s|ai)} 
where l/(|s|ai) = oo when ai = 

Proof. Take the first and last terms of the sum (|8.ip to get 

p(n) n — 2 

an> — r + aiCLn-i (8.2) 

for n > 2. In particular 

p{n) 



an > 



so Pa < P- Further if oi and at least one a^, j > 1 are nonzero then by (j8.2p a„ > for all n > j, 
since the p{n) are assumed nonnegative. In this case, then, we also have 



(n — 2)ai|s| 

and so pa < l/(|s|ai). The inequality pa < l/(|s|ai) also holds by convention if ai = 0. Finally 
suppose ai 7^ but a„ = for all n > 1. Then p(n) = for all n > 1, and, from (jS.ip for n = 2, 
s = 1. This is the case we have excluded. The result follows. □ 

For the lower bound on the radius we need a few preliminary results. First, some simple 
combinatorial facts. 



Lemma 8.4. 

for n,k £ Ij, n > k > 0. 
Proof. 



k - \kJ 



n\ nn — 1 n— k+l^nn n (n 



k I k k — 1 1 k k k 



□ 



1 /n\ 6'-J+i 
> 



Lemma 8.5. Given < 6 < 1 



n \3J J 

for 1 < j < On and n> 2. 

Proof. Fix n. Write j = An, < X < 0. Then using Lemma 

1 /n\ _ 1 / n \ n^"-i _ A^^^+i ^ 0-^+^ 
"nxij n\XnJ ~ (An)'*'" An ~ j 



□ 
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Second, we need to understand the behavior of ^ a„3;" at the radius of convergence. 
Lemma 8.6. 

A(x) < f(x)+x\s\(l + e)A'(ex)A(x) + (x^A^(e^x)) + PJx) 

2x ax \ / 

for all < 6 < 1/e, e > 0, and < x < pa, where Pe{x) is a polynomial in x with nonnegative 
coefficients. 

Proof. Take < 6* < 1/e and e > 0. 



{n) (, ,n . , . \ ^ / n\ ^ 
— + \^s\- - sign(s) J 2^ J a^an-j 



<^ + |.|(n + 2) E J a,a„_^ + kl^ E J 

l<j<en ^-^^ en<j<n-en ^"^ ^ 

by Lemma 18.51 

.p(n) ,n + 2 \s 



l<j<en en<j<n-en 

by Lemma 18.51 and since (x/n)^ is decreasing for < x < n/e 

Thus for n sufficiently large that (n + 2)/n < 1 + e the coefficients of A(x) are bounded above by 
the coefficients of 

fix) + x\s\{l + e)A'(6'x)A(x) + ( x'^A^O'^x) 

2x ax V 

Adding a polynomial to dominate the earlier coefficients of A(x) we get that the coefficients of 
A(x) are bounded above by the coefficients of 

fix) + x\s\{l + e)A'{ex)A{x) + '^-^^ (x^A^{9'^x)] +P,{x). 

2x ax \ ) 

Since all coefficients are nonnegative, for any < x < we have 

A(x) < /(x) + x|s|(l + e)A'(6'x)A(2;) + (x^A^{B\) \ 

2x ax V / 

□ 

Lemma 8.7. If pa < p and pa < l/(|s|ai) then A{pa) < oo. 
Proof. Consider Lemma 18.61 Choose 9 > and e > so that 
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which is possible since limg^Q A'{9x) = oi and pa < l/(|s|ai). Letting x ^ pa we see that 

hm A{x) < C + pa\s\{l + €)A'{9pa) hm A(x) 

where C is constant, since 9^ < 1, and pa < p- So 

(1 - pa\s\il + e)A'{9pa)) hm A(x) < C. 

But by I- pa\s\{l + e)A'{9pa) > 0, so A{pa) < oo. □ 

Lemma 8.8. If pa < p and pa < l/(|s|ai) then A(x) is unbounded on < x < pa- 
Proof. Take any e > 0. Then there exists an > such that for n > N 

p(n) 

n! ^ 



Define 



■r ^ pM ill I v^n— 1 

an if On > ^ + \s\CiCn-l + £ 2_^j=i CjCn-j 

+ |s|ciCn-i + CjCn-j otherwise (in particular when n > N) 



In particular ci = ai. Let C(x) = X]a,>i CnX*^ (which implicitly depends on e) have radius p^- Since 
o-n ^ Cn, Pa ^ Pe- Rewriting with generating series 

C(x) = f{x) + \s\aixC{x) + eC^{x) + Pe(a;) 

where -Pe(2;) is some polynomial. This equation can be solved by the quadratic formula. The 
discriminant is 

A, = (1 - \s\aixf - 4e{f{x) + P,(x)). 

Pe is the closest root to of A^. 

By construction, the coefficient of x" in P^ix) is bounded by a„. Suppose A{pa) < oo. Thus 

f{Pa)+Pt{Pa) < f{Pa)+-^{Pa)- By the nonuegativity of the Coefficients of / and then /(x) +Pe (2;) 
independently of e for for Q < x < pa- Thus 

lim Ae = (1 - \s\aipaf 

for < X < /9a. So 

1 1 

> Pa> Pe 



\s\ai \s\ai 

as e ^ which is a contradiction, giving that A(x) is unbounded on < x < p^. □ 

Proposition 8.9. pa > minjp, l/(|s|ai)}, where l/(|s|ai) = oo when ai = 0. 

Proof. Suppose on the contrary that Pa < P and pa < l/{\s\ai) then Lemmas 18.71 and [8 .81 contradict 
each other so this cannot be the case. □ 

Taking the two bounds together we get the final result 

Theorem 8.10. Assume Ylk>i-^^Pi^)/^^- radius p. Then X]^"7i,™/'^' converges with radius 
of convergence min{/7, 1/(571,1)}, where l/(|s|7i,i) = co i/71,1 = 0. 

Proof. Immediate from Lemmas 18.31 and 18. 9i □ 
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8.2 Systems 

Now suppose we have a system of Dyson-Schwinger equations as in ([37 



X^{x) = I - sign(.,) J2 ^'B^^-'^iX^Q') 



k>l i=0 

for r £ TZ with TZ a finite set and where 



X} '''' 



for all r G TZ. 

To attack the growth of the 7[ we will again assume that the series of primitives is Gevrey-1 
and that the Sr give nonnegative series. 

Assumption 8.11. Assume s,- > 1 or s,- < for each r £ TZ. Assume that 

k>l 

has radius < ^ oo, p'^'{k) > for k >1, and the /''(x) are not identically 0. 

We'll proceed by similar bounds to before. 

Definition 8.12. Let < = 7[^„/n! and A''{x) = X;„>i a'i^x''. 

Again the a[ are all nonnegative. 
Then 

< = ^ + J2^\^r\^-sign{sr))ala:Jfj ' + ^ X^d^^l^K-^ (l) ^ (^.4) 
i=i ^ ^ j&n i=i ^ ^ 

Proposition 8.13. For all r G TZ, the radius of convergence of A'^{x) is at most 

minL,,— 

interpreting the second possibility to be oo when ^j^n \s,\a{ = Q. 
Proof. Taking the last term in each sum of (j8.4p we have 

P' {n) I j \ n — 2 



n! \ ^ ■ - ■ - / n 



Let 6^ be the series defined by b^ = Ui and equality in the above recursion. Then argue as in the 
single equation case, Proposition 18. 3^ to get that the radius of A^(x) is at most 



mm < pr 



□ 
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Proposition 8.14. The radius of convergence 0/ ^^.^^ A''(x) is at least 

/ 1 1 

mm < Pr , r > 

interpreting the second possibility to be 00 when '}2j<^n 

Proof. The overall structure of the argument is as in the single equation case. 
The equivalent of Lemma 18.61 for this case follows from 

n-l 





^ p'{n) ^ 


n + 2 ^ 1 
n ^ 




rg7^ 










~ ^ n\ 


n + 2 , 

/ -s? 

n ^ ' 


\a{ Y fln-i 












+ max( 
i 


n-2 , 








1=2 ^ 


^ \re7^ 











+ max(|.,|)(n + 2) (E<-)(E<) 

for 9 as in Lemma r8.6l with X^^^t^ A''(3;) in place of A(x), where A''(3;) = ^ a''(n)x". Then continue 
the argument as in Lemma 18.61 with Ylrenf^i-^) ™ place of f{x) and maxj(sj) in place of s, and 
using the second term to get the correct linear part. 

For the argument as in Lemma 18.81 Take any e > then there exists an > such that for 
n > N we get 



Taking C''(a;) to be the series whose coefficients satisfy the above recursion with equality in the 
cases when this gives a result > and equal to aj^ otherwise we get 

\jeTi J j&n 

where PJ is a polynomial. 

Summing over r we get a recursive equation for ^^^^7^ C^(x) of the same form as in the single 
equation case. Note that since each is a series with nonnegative coefficients there can be no 
cancellation of singularities and hence the radius of convergence of each C is at least that of the 
sum. Thus by the analysis of the single equation case we get a lower bound on the radius of 

A''(2;) of m.msen{Ps, '^/T. jan \sj\ai]- □ 
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Proposition 8.15. Each A^{x), s £ TZ, has the same radius of convergence. 

Proof. Suppose the radius of A^{x) was strictly greater tlian tliat of A''(x). Tiien we can find 
(3 > 5 > such that 

a; > /5" > 5" > 
for n sufficiently large. Pick a k > 1 such that > 0. Then 

^ " - n---(n-/t + l) "-'^ ^ n---{n-k + l)^ 

so 

> 



|sr|a| n - ■ ■ {n — k + 1) 

which is false for n sufficiently large, giving a contradiction. □ 

Theorem 8.16. For all r gTZ, X]^'^7in/"'' converges with radius 



rev?. ^ — ' ' 



where the second possibility is interpreted as oo when Ylji^^n \ 

Proof. Take s £ TZ such that ps is minimal. 

Since we are working with nonnegative series the radius of A'^(a;) is at least that of ^j-eiz 
Hence by Lemmas 18.131 and 18.141 A^(x) has radius exactly 



min{/)^,l/ V|sj|a{}. 

re7<. ^ — ' 

Thus by Lemma [8. 151 all the '^af^x^ have the same radius minrg7^{/9r, 1/ X^^^t^ |sj|a-{} □ 

8.3 Possibly negative systems 

Let us relax the restriction that p'^{n) > 0. It is now difficult to make general statements concerning 
the radius of convergence of the A'"(x). For example consider the system 

If N n-1 , . _i 71-1 . s _i 



n-l / \ "^^^ / \ 



so si = 2 and S2 = —1- Suppose also that 



p\2) = -Aiajf 

1 2 

= ai 

p\n) = -2{n-iy.alai_^ 
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Then = and inductively = for n > 2 so the system degenerates to 



n— 1 



ai=^ + ^(2j-l)aX 



n-j 



n\ ^ n — 1 I 1 



n I a| if n = 1 

«n = S 

I otherwise 

We still have a free choice of p^{n), and hence control of the radius of the series. On the other 
hand the series trivially has infinite radius of convergence. 

Generally, finding a lower bound on the radii of the solution series, remains approachable by 
the preceding methods while control of the radii from above is no longer apparent. 
Precisely, 

Theorem 8.17. The radius of convergence of ^^^^x^'jl^^/nl is at least 

1 



mm < 



Pr, 



where the second possibility is interpreted as oo when ^j^-ji i = 0- 
Proof, for any e > 

n-2 



\aj < ; h 



jen 

n-2 



4-i| + Yl Kl'^'-I* ~ sign(s^))||a- ||a;_ 



1=1 
-1 



n-l 



So, for a lower bound on the radius we may proceed as in the nonnegative case using the absolute 
value of the coefficients. □ 
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Chapter 9 

The second recursion as a differential 
equation 

In this final chapter let us consider the second recursion derived in Chapter [7] as a differential 
equation rather than as a recursive equation. That is, in the system case 

7[(x) = Prix) - sign(s,)7[(x)2 + \s,\-fi{x)^ xd^Yiix) (9.1) 

as r runs over TZ, the residues of the theory. While in the single equation case 

m7i(x) = P{x) — sign(s)7i(2;)^ + \s\ji{x)xdxji{x) (9-2) 

The parameter m was added to keep the QED example in the most natural form, however it is not 
interesting since since we can remove it by the transformation 71 (x) 1— > 77171 (x), P{x) 1— > m?'P{x). 

No non-trivial results will be proved in this chapter, we will simply discuss some features of 
some important examples. More substantial results will appear in |31j . 

As a consequence of the renormalization group origin of the first recursion discussed in section 
14.11 the /3-function for the system shows up as the coefficient of (7i)'(a;), namely 

I3{x) = \sj\l{{x) 

in the system case and 

j3{x) = x|s|7i(x) 

in the single equation case. Consequently this differential equation is well suited to improving our 
understanding of the /3-function. 

In particular in the single equation case, we see immediately from (|9.2|) that any zeroes of (3{x) 
must occur either where P{x) = or where 7[(x) is infinite. The second of these possibilities does 
not turn out to be physically reasonable as we will discuss in more detail below. The system case 
is not quite so simple. Assume (3{x) = 0. If we rule out infinite (7i)'(a;), then we can only conclude 
that for each r £ TZ 

Yi{x) + sign(s^)7];(x)2 - P,.{x) = 0. 

In order to extract further information in both the single equation and the system case we will 
proceed to examine plots of the vector field of (7[)'(x), first in some toy single equation cases, 
second in the case of QED reduced to one equation, and finally in the 2 equation example of cp'^. 
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S = 2 




Figure 9.1: The vector field of 7i(x) with s = 2, m = 1, and P{x) 



X. 



9.1 Toys 

First let us consider a family of examples which are simpler than those which occur in full quantum 
field theories, namely the family where m = 1 and P{x) = x. 



9.1.1 The case s = 2 

If we set s = 2 we get the situation explored in [S] which describes the piece of massless Yukawa 
theory consisting of nestings and chainings of the one loop fermion self energy into itself as discussed 
in Example I3.5[ The second recursion viewed as a differential equation is 

71 (x) = X - 7i(j;)(l - 2xdx)-fi{x). 

Broadhurst and Kreimer [5] solved this Dyson-Schwinger equation by clever rearranging and 
recognizing the resulting asymptotic expansion. The solution, written in a slightly different form, 
is given implicitly by 

with integration constant C. 

We can proceed to look at the vector field of 7i(x), see Figure EH 

We are primarily interested in the behavior in the first quadrant. Of particular interest are 
possible zeros of solutions since, in this simple single equation situation, x'ji{x) = f3{x) where (3{x) 
is the /3-function of the system. 
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Figure 9.2: Solutions which die in finite time along with the curve where 7i(x) = 0. 
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Figure 9.3: The four loop approximation near the origin. 
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From the figures we notice a family of solutions which come down to hit the x axis with vertical 
tangent. These solutions have no real continuation past this point. These solutions are consequently 
unphysical. It is not clear from the figure whether all solutions have this behavior. One of the 
major goals of [31] is to find conditions guaranteeing the existence of a separatrix. 

Viewing the vector field near the origin can be quite misleading, since it appears to have both 
types of behavior simply because all the solutions have the same asymptotic expansion at the origin. 
Additionally this implies that the apparent, but potentially false, separatrix is well matched by the 
first four terms of the asymptotic expansion as illustrated in Figure 19.31 Of course given that we 
have a recursive equation and an implicit solution we can easily calculate the asymptotic series out 
to hundreds of terms [5], and the use of a four loop approximation is merely meant to be illustrative. 

Another simple observation is that we can derive the equation for the curve where the solutions 
are horizontal by solving for 7j(x) 

// X ii{x)+iKx) -X 

' 2x71 (x) 

and then solving the numerator to get the curve 

-1 + yJl + Ax 

y= — 2 — 

illustrated in Figure [9T2l 
9.1.2 Other cases 

Let us return to general s while maintaining the assumption m = 1, P{x) = x. 

The case s = is degenerate, giving the algebraic equation 71 (x) = x — 71 (x)^ with solutions 

= ^ • 

From now on we will assume s ^ 0. 

We can obtain implicit solutions for a few other isolated values of s using Maple 

s = 1 :7i (x) = X + xW ( C exp ( 

s=^:A{X)- xi/32i/3^' (X) =C^B {X) - xi/321/35' (X)) where X = ^jj^- 

s = 3 :(7i(x) + 1)^ (X) - 2^/3^' (X) = C ((71 (x) + l)B (X) - 2^/3^' (X) 

(1 +7i(x))2 + 2x 



where X 



24/33.2/3 



where A is the Airy Ai function, B the Airy Bi function and W the Lambert W function. 

Qualitatively the vector fields are rather similar, see Figure 19.41 The same qualitative picture 
also remains for values of s > where we do not have exact solutions. For s < the picture is 
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s = 1 



s = 1.5 
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(b) s = 3/2 
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(c) s = 2 (d) s = 3 

Figure 9.4: The vector field of ^'\{x) with m = 1 and P{x) = x, showing the dependence on s > 0. 
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Figure 9.5: The case P{x) = x and s = —2. A typical example with s < 0. 
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Figure 9.6: The case s = 1 compared to the curve 71 (x) = x 
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somewhat different, see Figure 19.51 but we still see solutions which die and can still ask whether 
there are solutions which exist for all x > 

In the case s = 1, 71 (x) = x is manifestly a solution, so there are solutions which exist for all x 
for some values of s. 71 (x) = x is illustrated in Figure WM 

Note also that we can, as before, calculate the curve where solutions are flat for general s, and 
it depends only on the sign of s since 

1, . 7i(x) + sign(s)7^(x) -x 

lA^) = n — n 

|s|x7i(xj 

giving the curve 

— 1 + -v/l + sign(s)4x 

y= 2 • 



9.2 QED as a single equation 

In this section we are interested in the case where m = 2, and s = 1 in (j9.2p . In view of the Ward 
identities and the work of Johnson, Baker, and Willey [19] the QED system can be reduced by a 
suitable choice of gauge to the single equation with those values of m and s describing the photon 
propagator. 

The first question is how to choose P{x). To 2 loops 

2 

To 4 loops we need to correct the primitives in view of the reductions of the previous chapters. 
Values are from [16]. 

2 

P[x) = I + + (-0.0312 + 0.06037)x3 + (-0.6755 + 0.05074)x^ 

In the first of these cases little has changed from the simple examples of the previous sections. 
At 4 loops, however, P(0.992 . . .) = which causes substantial changes to the overall picture, see 
Figure 19.71 

This zero in P{x) is expected to be spurious, due only to taking the 4 loop approximation out 
beyond where it is valid, and the qualitative behavior of the solutions looks much more familiar if 
we restrict our attention to < x < 0.992 . . ., see Figure [9181 

Note that if P{x) > for x > then by the same analysis as in the P(x) = x case we can 
determine the curve where the solutions are flat. The curve is 



-1 + V^l + 4P(x) 

y= 2 • 

The first four loops of perturbation theory give a good approximation to reality, and also as 
expected match the apparent separatrix for small values of x, which is illustrated quite strikingly 
in Figure \9M 
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QED to 4 loops 
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(a) P{x) taken to 2 loops. 



(b) P{x) taken to 4 loops. 



Figure 9.7: The vector field of 7((x) for QED with different choices for P{x) 
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Figure 9.8: The region between x = and x = 1 in the vector field of 7((x) for QED with P{x) 
taken to 4 loops. 
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(a) Close to the origin, with x = 1/137 . . . marked. 



(b) The limits of the validity of the four loop approx- 
imation. 



Figure 9.9: The four loop approximation to 71 (x) for QED. 



9.3 04 

Let us now consider 0^ as an example which legitimately leads to a system of equations, but for 
which it remains possible to create illustrations, and perhaps even to analyze. Taking advantage 
of the graphical similarity between the vertex and propagator in cp^ and the symbols + and — 
respectively we will write the specialization of ()9.ip for (p'^ as the system 

^+{x) = +7+(x)2 + (7+(x) + 27f (x))xa^7+(x) 

ry-(x) = p-{x) - 7f (x)^ + (7^(x) + 27f (x))xaa;7f (x) 

The values of 7^*" and 7^" up to order x^ can be obtained from [20j and hence so can those of 
and P~ . Close to the origin we see a distinguished solution, see Figure [9T0l As in subsection 19. 1.1^ 
this may not indicate a solution which exists for all x, but we hope that this solution is physical. 

There are many tantalizing features appearing in these examples which will hopefully be the 
genesis for future work linking to different fields. The equations derived in Chapter [7] seem consid- 
erably more tractable than the original Dyson-Schwinger equations when viewed either as recursive 
equations or as differential equations. They have already led to physically interesting results as in 
Chapter [8] and hold much promise for the future. 
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